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CHAPTER I . SUMMARY 


It seems possible that, in supersonic flight, unconventional arrange- 
ments of wings and bodies may offer advantages in the form of drag reduc- 
tion. It is the purpose of this report to consider the methods for deter- 
mining the pressure drag for such unconventional configurations, and to 
consider a few of the possibilities for drag reduction in highly idealized 
aircraft . 

The idealized aircraft are defined by distributions of lift and 
volume in three-dimensional space, and Hayes' method of drag evaluation, 
which is well adapted to such problems, is the fundamental tool employed. 
Other methods of drag evaluation are considered also wherever they appear 
to offer simplifications. 

The basic singularities such as sources, dipoles, lifting elements 
and volume elements are discussed, and sane of the useful inter-relations 
between these elements are presented. Hayes' method of drag evaluation 
is derived in detail starting with the general momentum theorem. 

In going from planar systems to spatial systems certain new problems 
arise. For example, interference between lift and thickness distributions 
generally appears, and such effects are used to explain the difference 
between the non- zero wave drag of Sears -Haack bodies and the zero wave 
drag of Ferrari 's ring wing plus central body. 

Another new feature of the spatial systems is that optimum configu- 
rations generally are not unique, there being an infinite family of lift 
or thickness distributions producing the same minimum drag. However it 
is shown that all members of an optimum family produce the same flow 
field in a certain region external to the singularity distribution. 

Other results of this study indicate that certain spatial distri- 
butions may produce materially less wave drag and vortex drag than com- 
parable planar systems . It is not at all certain that such advantages 
can be realized in practical aircraft designs, but further investigation 
seems to be warranted. 


* 


Unedited by the NACA (the Committee takes no responsibility for the 
correctness of the author's statements). 
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CHAPTER II. INTRODUCTION 


The primary purpose of this report is to consider the problems 
involved in exploring a broader class of aircraft configurations than 
is ordinarily studied. for supersonic flight. It is necessary to deter- 
mine whether any unconventional arrangements of wings and bodies offer 
sufficient aerodynamic advantages in the form of drag reduction to merit 
more detailed study. As a first step in this direction attention is 
directed to optimum configurations, even though they are highly idealized 
in form and do not necessarily represent practical aircraft. 

In the preliminary exploration of such configurations it is not 
necessary to know their detailed shapes. It is stiff icient to define the 
aircraft as a distribution of lift and volume in space, without knowing 
the camber and twist of the wing surfaces supporting the lift distri- 
bution, and knowing only approximately the shapes of -the bodies con- 
taining the volume. 

Hayes' method of drag evaluation is well adapted to this type of 
analysis and is one of the primary tools used. However other methods 
and points of view are employed wherever they appear to offer further 
understanding of the problems. ... — 

The properties of sources, dipoles, etc., are reviewed, and a sin- 
gularity corresponding to an element of volume is introduced. Some 
useful relations between three-dimensional distributions of different 
types of singularities are developed and later applied. Also Hayes’ 
method for drag evaluation is developed in detail. 

l 

Since this report is exploratory in nature the investigations made 
are frequently incomplete and somewhat isolated from each other. Some 
of the material of Ref. 2 and most of the material of Ref. 3 are included 
in this report for convenience. The latter has also been published in 
The Aeronautical Quarterly, May 1955, under the title, "The Drag of 
Non-Planar Thickness Distributions in Supersonic Plow." Permission to 
reproduce this material has been granted by The Royal Aeronautical 
Society. 
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CHAPTER III. SINGULARITIES UT TT.T7.E n IN THE "LINEARIZED" 


DESCRIPTION OF THE FLOW ABOUT AIRCRAFT 


A. BASIC SINGULARITIES 


The Source 


For incompressible, non-viscous fluids the equation governing the 
flow is the Laplace equation. 


ii. + ii + ii. , o 

cix 2 dy 2 dz 2 


(3a-l) 


where (f> is the perturbation velocity potential. A basic solution, ■which 
exhibits spherical symmetry, is the source. 


0 S = t ===== === ^ L============ (3a-2) 

1 M(x - |) 2 + (y - n ) 2 + (z - £) 2 


■i 


This solution can "be interpreted an representing the emanation of unit 
volume of fluid per unit of time from the point r|, Because of 
the linearity of Eq. (3a-l), other solutions of it can he huilt up by 
a superposition of sources through the use of certain limiting proce- 
dures; such resulting solutions are the horseshoe vortex, doublet, line 
vortex, etc. Much is known about these solutions and with them the 
flow over wings and bodies can be described mathematically. 

In supersonic flow the governing differential equation is the 
linearized potential equation. 


-0 


cbc^ Sy2 dz^ 


(3a-3) 


where x is the coordinate in the stream direction and p -Hvr - 1. Equa- 
tion (3s-3) can also be considered as the two-dimensional wave equation 
where the x coordinate is thought of as the "time" variable. 

If tbe y and z coordinates in the Laplace equation (Eq. 3a-l) are 
multiplied by ip, then that equation is transformed into the wave equa- 
tion; a similar transformation of the source potential from Eq. (3a-2) 
results in 
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(3a-4) 


which can easily he shown to he a solution of Eq. (5a-3). Equation (ja-4) 
is real inside the forward and rear Mach cones, (x - | f > 0 2 |?y - „) 2 + 

(z - £ ) 2 J , and imaginary elsewhere; however, due to the nature of super- 

sonic flow only the solution in the rear Mach cone is used to represent 
a source. Since half of the real solution is discarded, the constant 



associated with the incompressible 0 must he doubled to represent a unit 
supersonic source. Thus the supersonic source at g, tj, £ has the potential 


r 


-1 

■*-5 ^pf(y-n) 2 + (a- t) 2 

(3a-5) 

2rtf(x -|) 2 - p 2 [(y -n) 2 + (z- £) 2 ] 

0 

Elsewhere 


where the x axis is in the free stream direction. It can he shown that 
Eq. (3a- 5) represents unit volume flow from the point §, rj, £; however, 
care must he taken in the proof because of the singularities on the Mach 

cone (cf. Ref. 4). In the proof given by Robinson he made use of the 
concept of the finite part of an infinite integral, an idea originally 
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(c 5 ) 

introduced "by Hadamard v>; . As in incompressible flow, other solutions 
of Eq.. (3a— 3) can be built up by superposition of the basic source solu- 
tions; some of the solutions can also be obtained, as waB the source, 
by analogy with the incompressible solutions. 

Before going on to other solutions let us examine the supersonic 
source in more detail. Since the velocities are infinite on the Mach 
cone from a finite source, care must be taken in using such sources to 
describe real flows. It is instructive to examine the isolated source 
in terms of the limit of a finite line of sources in the free stream 
direction as the length tends to zero while the total strength remains 
constant. Under the assumptions of slender body theory, if the line of 
sources extends from x = 0 to x = x Q with strength Kx, it represents a 

cone of semi -vertex angle K/2U with a semi-infinite cylindrical after- 
body (Pig. 3a-la) . The velocities are constant along conical surfaces 



Fig. 3a— la: Cone-cylinder and Fig. 3a— lb: Perturbation flow 
source distribution repre- lines in x-z plane for super- 

senting it sonic source, M = ^2 
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from the origin; hut on the Mach cone from x = x Q the velocities become 
infinite due to the discontinuity in source strength. The total inte- 
grated source strength C is equal to l/2 Kxq^. If xo is allowed to 
approach zero while C remains constant then in the limit a concentrated 
source of strength C is obtained. The flow pattern in the xz plane for 
the source at M = {2 is shown in Fig. 3a.-lb_. (See also Ref. 6.) For a 
source of finite strength the velocities are infinite on the Mach cone. 


The Three-Dimensional Doublet 

The three-dimensional doublet (or dipole) is a second basic solution 
of the wave equation; it is obtained by allowing a source and sink of 
equal strength to approach one another while the product of source strength 
and distance between source and sink remains constant (and equal to unity 
for a unit doublet). The axis of the doublet is defined here as the 
vector extending from the cehter of the sink to the center of the source; 
positive values are taken to be those along the positive directions of 
the coordinate system. For a doublet with its axis vertical, the above 
method of derivation Is equivalent to taking the negative partial deriva- 
tive of the source potential in the z direction; that is. 



x > |Br 


(3a -6) 


where r 2 = y 2 + z 2 . Equation (3a -6) represents, a positive doublet at 
the origin, l.e., one with the source above the sink. 


The Horseshoe Vortex - 

In supersonic theory, as well as in subsonic, the_flow around a 
wing of finite span can be described by certain solutions of the wave 
equation called horseshoe vortices. In the subsonic case this singularity 
Is derived by Integrating in the streamwise direction a semi-infinite line 
of negative doublets with axes vertical. The supersonic horseshoe vortex 
can be derived in the same way as the subsonic one If only the finite 

part of the integral, as defined by Eadamard^^ , is taken as the solution. 
This solution can also be obtained without the use of the Hadamard finite 
part if a streamwise line of sources is dif f erentiated . in the vertical 
direction; thus, 
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x :> 0r 


(5a-7) 


The flow pattern for a horseshoe vortex in planes normal to the 
free stream axis is shown in Fig. 3&-2. Far behind the bound vortex, 


\ 


4 



Fig. Ja- 2 : Flow pattern for supersonic horseshoe vortex 
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the flov near the x axis is similar to the flow far downstream, around 
a subsonic horseshoe vortex and it is this part which gives rise to the 
"vortex" drag. The drag associated with. the flow near the Mach cone 
is called "wave" drag. Equation (?a-7) represents a supersonic horseshoe 
vortex of unit strength, i.e., unit circulation around the hound vortex. 
Since a force pUr is associated with a hound vortex of strength r, we 
shall, for convenience, discuss unit lifting elements which have as their 
velocity potentials 


*x 


xz 

2rtpUr 2 l/x 2 - 


x > ^r 


(3a-8) 


Similarly, the potential for a unit side force element is 



x > pr 


(5a -9) 


The force associated with Eq. (3a- 9) is directed in the positive y direc- 
tion; a force in' any direction normal to the flow direction may he repre- 
sented hy a combination of lift and side force elements . 

In the light of the discussion of the horseshoe vortex, the three- 
dimensional doublet (Eg.. 3a-6) may he given added significance as a lift 
transfer element or element of moment. That is, the doublet potential 
can he formed hy subtracting the potential for a horseshoe vortex at x = Ax 
from one at x = 0 (Fig. 3 a ~3) and applying the proper limiting processes 
(equivalent to differentiating the horseshoe vortex potential.); in tills 
process the trailing vortices from the negative or rear element are can- 
celed out hy those from the positive one, and the remaining part forms 
the doublet or . lift transfer element . - 

Z Z z 




Fig. 3 a -3: Formation of doublet or lift transfer element from 

horseshoe vortices 
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The horseshoe vortex consists of a bound vortex of infinitesimal 
length plus two free vortices trailing back to infinity. Since the 
vortex drag and the lift associated with a finite wing can be evaluated 
by considering the flow velocities far behind the wing, it is useful to 
consider the trailing vortices as they appear in the Trefftz plane far 
downstream from the bound vortex. The Trefftz plane flow represents a 
two-dimensional doublet or dipole and its potential is obtained by 
letting x — ><*> in Eg. (3a-7). Thus 



z 

2rt(y 2 + z 2 ) 


(3a-10) 


It should be noted that the potential for this doublet is independent 
of Mach number, and thus the vortex drag calculations for a given lift 
distribution are the same for supersonic and incompressible flows. The 
flow pattern about the doublet will be similar to the planar flow inside 
the dashed circle in Fig. 3&-2. 


The Volume Element 


Another useful solution is the doublet with its axis in the stream- 
wise direction; it has as a potential 


= - 



; oc 


x > (3r 


(3a-ll) 


Equation (3a-ll) can be shown to represent the potential for a unit of 
volume equal to l/u (see Chapter VTI) at the origin. A distribution of 
volume elements along the x axis with strength f(x) has as a potential 



f(i)0 v <34 


f(i) 

x-(Br ox-0 r 

f'(6)d| 

2 *f(x - |) 2 - P 2r 2 

0 2rt ^ 0 

f(x - |) 2 - 


(3a-12) 


The first term in Eq. (3a-12) is infinite; but if only the finite part 
of the integral is considered ^as defined by Hadamard^^, then Eq. (3a- 12) 
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represents the potential for a source distribution of strength f '(§)• 
Thus a body can he built up from, a series of volume elements as well as 
from a series of sources and sinks. 


The Closed Vortex Line 

Equation (3a-ll) can be considered not only as a volume element but 

also as a closed vortex line of circulation strength l/g^ in the yz plane 
(Fig. 3a-4a). The line carries a constant intensity of forces directed 
inward so that the total vector force is zero. The negative of Eg. (3e— 11) 
would represent an element with the forces directed outward from it. The 
potential for the closed vortex line can also be obtained by applying the 
standard limiting process to an element composed of two pairs of horseshoe 

vortices of strength l/|3^, one with its axis in the negative z and the 
other in the negative y direction (Fig. 3a— 4b); when added together the 
trailing vortices cancel leaving the closed vortex line . 



Fig. 3a-4: Formation of closed vortex line from horseshoe vortices 


Two-Dimensional Singularities r; 

In subsonic flow two-dimensional sources, obtained by integrating 
a line of three-dimensional sources in the lateral direction, have proven 
useful in many problems; so also has the infinite bound vortex obtained 
by a lateral integration of horseshoe vortices. The same types of solu- 
tions can be derived for supersonic flow and these provide more insight 
into the nature of the supersonic solutions. The two-dimensional source 
potential is 
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02S = J' 


w- 


p2 z 2 


;2-p2 z 2 


0S dT l - - h x ^ 


2P 


pz 


(5a- 13) 


All of the disturbance created by the two-dimensional source .is 
concentrated on the Mach planes from, it, thus creating a potential jump 
across these planes. The two-dimensional vortex potential is 


02V - 


|/x^-p2 z 2 


0HSV <3-n = 


+ i 

z > 0 

1 

" 2 

z < 0 


x> |pz| (3a-l4) 


Again all of the disturbance is concentrated on the Mach planes. There 
is a potential jump across the Mach planes and also across the z = 0 
plane, the latter due to the discontinuity in the past history of the 
fluid particles above and below the plane. 
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B. SOME EQUIVALENT SINGULARITY DISTRIBUTIONS 


Statement of the Problem 

The first section of this chapter reviewed the basic singularities 
which represent elements of lift, side force and volume in linearized 
supersonic flow theory. It was noted that these singularities may all 
be obtained from the source singularity with the aid of the simple proc- 
esses of integration and differentiation. The fact that the basic 
singularities are so related will be shown to imply that certain dis- 
tributions of singularities are equivalent, i.e. they produce the same 
flow field, at least outside of a finite region. In._the present section 
an equivalence theorem will be proved regarding constant strength dis- 
tributions of sources, lifting and side force elements and vortex sheets. 
Such a theorem will later prove useful in the study of interference 
between distributions (Ch. IX B,C,F). Note that if the distribution A 
is part of a larger distribution (A,B) and if A is replaced by an equiva- 
lent distribution A^_ then the drag of ^Aj_,b) is the same as that of (A,B). 

This follows from the fact that the substitution of A-j_ for A does not 
change the flow field at infinity (Ch. IV). 

The distributions to be studied will be located on a cubic shell 
which has two faces perpendicular to the free stream direction. One 
face of the cube will be covered by sources of constant strength and 
the opposite face by sinks of constant strength. The remaining four . 
faces will be encircled by vortex lines of constant strength. Two cases 
may then be distinguished: (a) The source distribution is on a face 

parallel to the free stream; (B) The source distribution is on a face 
normal to the free stream. These two cases are illustrated in Fig. Jb-l. 
The source, sink and vortex distributions are uniform and of constant 
intensity as indicated. The vortex lines are continuous around the cube 
with the circulation directed so as to induce in the interior of the 
cube downwash velocities in case A and upstream velocities in case B. 



Fig. 5b-l 
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We shall now prove the following theorem. 


Theorem 

In “both, cases A and B the perturbation velocities are zero every- 
where outside the cubic shell. Inside the shell the downwash is constant 
in case A (w = -k) and the pressure is constant in case B (u = -k) . 

This theorem implies in particular that the source sink distribution, 
say In case A, is equivalent to the negative of the vortex-line* distribu- 
tion in case A, in the sense that the associated flow fields are Identical 
outside the cube. Note that in case A the vortex distribution on the 
front and rear faces gives rise to a lifting force, whereas the vorticity 
on the side faces produces no force. In case B the vorticity on the top 
and bottom produces lift and the vorticity on the side faces produces 
side force. 

The theorem will first be proved by a geometrical argument and then 
an alternative proof by analytical methods will be outlined. 


Geometrical Proof of Theorem 


Consider first cause A. We shall construct a geometrical configu- 
ration which corresponds to the distribution of singularities indicated 
in Fig. 3b-l. This construction will proceed in several steps by succes- 
sively cutting down configurations of infinite extent. The vortex dis- 
tribution on the front face is equivalent to a distribution of lifting 
elements of constant strength. 

To begin with we shall assume the whole infinite plane containing 
the front face to be covered by lifting elements. This may be physically 
realized by a cascade of doubly infinite (two-dimensional) wings of con- 
stant angle of attack a and such that the vertical dista nce betw een two 

neighboring wings is equal to the wing chord divided by / - 1, 

(Fig. 3b-2) . In the limiting case of zero chord length the plane x = -x 
is then covered by vortex lines with the circulation (of -strength k) 

oriented as in Fig. 3b-l* The value of the constant k is then k = 2aU. 

o 

The lift per unit area in the plane x = -3? is then 2apU . Since the 
wings are spaced so as not to interfere with each other hut still influence 
every point downstream of the cascade, the flow field at any point P with 
x > -x may he described as follows (Fig. 5b-2) . The point P receives a 
unit of downwash (-w = cdJ) frcoi the wings A and B each. It also receives 

a positive unit of pressure f-u = aU /]/ m^ - l) from A and a negative unit 
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Fig. Ja - 2 


of pressure from B. The net pressure preferred to p^j received at P is 
then zero and the net downwash is -w = 2aU = k. 

The cascade may now he terminated from above by a wedge of opening 
angle 2a located in the plane z = z" with its exterior surface parallel 
to the free stream direction (Fig. 3b-3) . Actually this wedge corre- 
sponds to a source distribution of constant source strength k = 2aU. 

If the cascade is removed for z > z the flow field is zero there since 
the wedge isolates this region from the rest of the cascade and since 
the exterior surface of the wedge is at zero angle .of attack. For z < z 
the flow field is unaffected by the introduction of the wedge. To see 
this consider a point P with z < z (Fig. 3^-3) • The wing at B acts as 
before to produce a downwash of -aU at P.’ Only the point C on the wedge 

VC 

' 

>P 1 ~ 

/ - 

/ 

B 


Fig. Jo-J 

affects the point P and this point C is already in the downwash field -aU 
of wing Ag. Thus the wedge turns the flow downward .only by an angle a 
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so that the total downwash at P is again -2aD. Conditions at P are the 
same as in the infinite cascade. 

Similarly the cascade may be terminated from below at z = -z by 
placing a wedge there of opening angle -2a. This corresponds to a dis- 
tribution of sinks of strength -k. 

The cascade may then be cut down to finite width by placing planes 
of zero thickness parallel to the plane at y = ty and removing the part 
of the wings for |y| > y. Since no sidewash is present the flow field 
is undisturbed by the introduction of these planes . Thus for x > -x, 

|y| < y, |z| < z the downwash is -w = 2aU = k and the pressure is zero. 
Outside this region all perturbation quantities are zero. 

Finally one may restrict the flow field to the inside of a cube by 
taking the negative of the above configuration and placing it at x = x. 

Thus the resulting flow field has constant downwash and zero pres- 
sure inside the cube -x < x < x, -y < y < y , -z < z < z. Outside this 
cube the perturbation velocity is zero. Thus the front face is a cascade 
of lifting wings at an angle of attack a, which bends the flow down. The 
rear face is a cascade of wings of angle of attack -a which straighten 
the flow out again. The top and bottom faces consist of wedges whose 
inside surfaces follow the direction of the flow which has been bent by 
the cascade. These outside surfaces are parallel to the free stream. 
(Note that for the wedge of negative angle the "interior" top surface is 
directed downward at an angle 2a and the "exterior" bottom surface is 
parallel to the flow. ) Finally the side faces are planes that carry no 
forces. For each such plane the downwash is -w = 2oCJ on the inside and 
w = 0 on the outside. These planes are then surfaces of constant vortic- 
ity. However, the vorticity vector is parallel to the free stream and 
hence no force results. 

Thus a geometric configuration (using a wedge of negative opening 
angle) corresponding to case A has been constructed and the theorem has 
been proved for this case . 

The corresponding construction for case B will only be indicated. 

The source distribution on the front face is obtained by placing wedges 
there (Fig. 3b-4) of half -angle a = 0k/2U. 


M r 





z 




Z 


-z 


Fig. 5b-^ 
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At a point P then the downwash is zero and the pressure is given by 
-u = 2ecu/^M 2 - 1 = k. 

By inserting planes of zero thickness at z = ±z, y = ±y and removing 
the wedges outside these planes the infinite configuration is cut down 
to a configuration with a finite cross section. Outside these planes the 
flow is undisturbed. Inside these planes -u maintains its value _ 

2aXJ^M 2 - 1 = k. These planes are then pressure discontinuities and 
hence carry lift and side force respectively. They _ are also vortex 
sheets. 

Finally the configuration may be terminated by placing its negative 
at x = +x. 

A geometric configuration (again using wedges of negative opening 
angles) corresponding to case B has thus been constructed and the theorem 
has been proved for case B. 


Analytical Proof of Theorem (Outline) 

Case A. Source Distribution Face Parallel to Free Stream 

Consider a cube with sources of strength k on the top and -k on the 
bottom, and with lifting elements of strength pllk on the front face and 
-pUk on the front face and -pUk on the rear face. On the side faces of 
the cube there are no forces associated with the vortex lines parallel 
to the flow direction; these are the trailing vortex system of the ele- 
ments on the rear face. 

In computing the potential due to the singularities on the cubic 
shell, various regions of the flow field are considered separately. For 
the region ahead of the foremost Mach waves from the cube no disturbance 
is possible in supersonic flow. Behind the cube, if the forward Mach 
cone from a point includes all of the shell, the potential at that point 
may be found simply by integrating the total effect of the singularities 
covering the shell. The potential due to individual^ unit source elements, 
lifting elements, and side force elements are given in Eqs. (ja-5)> (ja— 8), 
and (j5a-9). The strengths of the distributions considered in this case 
are indicated in Fig. 3b-5* — 
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Fig. 5b-5 

The potential for the entire shell is then 
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(3b-l) 


This, after evaluating the integrals, equals zero. 
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A third region of the flow field contains points slightly behind 
and far to the side of the cube, where forward Mach cones from the points 

include part, but not all, of the cube. For this region, Hayes' method 
can be used to show that the potential again is zero. This method is 
described in Ch. IVC. It requires that the distance from the cube to any 
point P where the potential is to be computed must be large compared to 
the dimensions of the cube . In addition, P must lie near the Mach cones 
emanating from the singularities on the cube. P is. then a point at acme 
angle 0 (measured from the horizontal plane) on a distant cylindrical 
control surface surrounding the cubic shell. An "equivalent lineal dis- 
tribution" of singularities is formed by finding the singularity strength 
intercepted from the cube by a set of parallel plane's originating at 
angle 0 on the control cylinder and inclined at the Mach angle to the 
free stream direction. The singularities intercepted by a given Mach 
plane are lumped together at the intersection of the Mach plane and the 
axis of the cylinder, such that the total strength of the equivalent 
distribution is equal to the total strength of the original distribution. 
After determining the strength (h) of the equivalent lineal distribution 
which represents the cubic shell for a fixed 0, the effect of all those 
singularities which influence the flow field at P can be summed. Hayes 
writes the expression for h as - 

h = +f - g z sin 6 - gy cos 0 

where f is the source strength (per unit length), g z j$ the circulation 
strength (per unit length) of the lifting elements, and Sy/P the cireu- 
• lation strength of the side force elements . 

Figure Jb-6 indicates the notation to be used in describing the 
geometry of the Intersections of the Mach planes with planes containing 
the x,y, z axes. The Mach plane is inclined to the axis of the control 

cylinder at the Mach angle p = sin“-*-(l/M) and it is tangent to a croBB- 
section of the cylinder at angle 0. The trace of the Mach plane in a 
horizontal (x-y) plane is inclined to a normal to the flow direction at 
angle 8, where tan 8 = cot p cos 0 = p cos 0 . The trace of the Mach 
plane In a vertical (x-z) plane forms an angle a with a line parallel 
to the z-axis, where tan a - 0 sin 0. 


(3b-2) 
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With this brief description of Hayes' procedure in mind, an equiva- 
lent lineal distribution of singularities is now to be computed for the 
specific case of the cubic shell described previously. Figure Jb-Ja. 
shows the intersections of two parallel Mach planes with the shell; the 
Mach planes are assumed to be separated by an infinitesimal distance. 

The case illustrated shows only three faces of the cube intersected by 
the Mach planes since the procedure would be the same if four faces were 
affected. In order to better define the geometry and notation, Fig. 3b-7b 
shows the cubic shell as though it were cut along the corner edges and 
flattened out in one plane . 
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The net singularity strength cut out by these Mach p lane s must be 
lumped along a length dx of the axis. The total source strength is the 
product of the strength per unit area (k) and the area intercepted from 
the top surface of the cube by the Mach planes: - - 


dn^ = k 


X 1 

tan & | 


(The negative sign is inserted because 9 is in the second quadrant for 
the example shown^ but f is positive . ) The total lifting element strength 
is pUk multiplied by the area intercepted from the front face: 
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Again, a negative sign is inserted because cos 0 is negative while l 
should be positive. There are no forces on the side faces. In computing 
the strength of the equivalent lineal distribution from Eq. (3b-2) it 
must be remembered that gy, g z from that formula are circulation strengths 

multiplied by 0 ; i.e.. 


®z 


PI 

pU 



Then 


h = +f - g z sin 0 - gy cos 0 


-kxj_ dx 
P cos 0 


3 I -pUkX! dx ' 
P u \p 2 sin 0 cos 0 , 


sin 0 


= 0 


(3b-4) 


That is, the net singularity strength is zero. This is true for all 
angles 0 , and similar calculations show that it is also true for every 
station x along the cylinder axis. Therefore, the velocity potential 
is zero at aJ i distant points for which Hayes' method is applicable. 

There remains to find the velocity potential in the neighborhood 
of the shell. The cube may be subdivided into smaller cubic shells, 
each similar to the original. Singularities on interfaces of adjoining 
shells then cancel so the net singularity distribution is unchanged. 

Those shells which lie behind and outside the forward Mach cone from 
any point cannot influence the velocity potential at that point. It 
was shown earlier that those shells which lie completely inside the for- 
ward Mach cone from the point also do not influence the potential there. 
Therefore, only those shells lying along the forward Mach cone need be 
considered. However, these may be further subdivided into cubic shells 
of elementary proportions so that the distance from the point to any 
one of the shells is very large compared to the dimensions of that shell. 
Then the analysis based on Hayes' procedure shows that these shells do 
not contribute to the velocity potential at the point either. This indi- 
cates that the velocity potential is zero everywhere outside the cubic 
shell. 

To find the velocity perturbations inside the shell, again consider 
it divided into smaller shells. None of these except the one containing 
the point P can influence the potential at P according to the preceding 
analysis. Therefore, all of the small shells located more than a dis- 
tance e ahead of P can be removed without affecting the potential at P. 
The forward Mach cone from P then intersects only the front face of the 
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remaining part of the original cube, so that, effectively, P is aware 
only of an infinite distribution of lifting elements. Since this result 
is independent of the location of P inside the original cubic shell, the 
downvash inside the shell must be constant. 


Case B. Source Distribution Face Normal to the Free Stream 

Consider now a cube with lifting elements of strength pUk on the 
top face and -pUk on the bottom face, with side force elements of 
strength pUk, -pUk on the side faces, and with sources of strength 
on the front face and -(3^k on the rear face. 


SOURCES 


L /FT/NG EL EMENTS 

'-SINKS 


(-XrJr?) 



SIDE FORCE 
EL E MEN Tsfp U-k ) 


Fig. 3b-8 


First, the potential ahead of the foremost Mach waves of the cube 
is, of course, zero. At a downstream point whose forward Mach cone 
includes ai 1 of the sources and lifting elements the__potential is 
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Carrying out the integration, it is found that 

(f = -2kx for -y < y < y and -z < z < z 
= 0 Elsewhere 
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By Hayes' procedure; when forward Mach, cones from distant points 
include only part of the singularities; the potential at those points 
is the same as would be contributed by a lineal distribution whose 
strength, h, can be computed in the manner described previously. For 
Mach planes intersecting the cube in the same location illustrated for 
another case in Fig. 3b-7; one finds that 


x-i dx X-, dx "" x-i dx 

f = _ k — r -= T -; S = -pUk -4——, l = -pUk 


sin 0 cos 0 


3 sin 0 


3 cos 0 


(3 b-6) 


and so 


h = -k 


X 1 


sin 8 cos 8 


cos 0 

pU\ (3 sin 0/ 


L 


x, dx \ 

— 1-pUk — sin 0=0 

pU \ 3 cos 0/ 


(3 b-7) 


Thus; the potential due to the cubic shell is zero at all distant points 
of the flow field which lie near the Mach cone of the shell. 

In the neighborhood of the cube, the same arguments used for the 
first cubic shell show that the perturbation velocities for this case are 
zero there also. Therefore, the perturbation velocities are proved to 
be zero in every region of the flow field external to the cubic shell. 

To find the potential at a point inside the shell, the shell is sub- 
divided as before into smaller shells, each similar to the original. The 
analysis just completed shows that the velocity perturbations at P cannot 
be influenced by any of these shells except the one containing P. There- 
fore, all of the small shells located more ;than a distance e ahead of P 
can be removed. The net singularity strength intersected by the forward 
Mach cone from P then includes only source ffTon the front face of the _ 
remaining group of cubes. Effectively, then, conditions at P are the 
same as behind an infinite distribution of sources of constant intensity. 
This result Is independent of the location of P inside the cubic shell, 
so the pressure must be constant inside the shell and the potential is 
of the form = cx . 
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CHAPTER IV. THE EVALUATION OF DRAG 


A. THE "CLOSE" AND THE "DISTANT" VIEWPOINTS 


The non-viscous drag for a wing and body moving at supersonic speeds 

may be obtained from two different points of view^), using linearized 
theory. First, the drag can be evaluated by integrating the local pres- 
sure times frontal area over the wing and body surfaces. Second, the 
drag can be evaluated from momentum or energy considerations involving 
the flow field at a great distance from the aircraft. These two pro- 
cedures are actually variations of the same basic method. 

In the latter case part of the drag due to lift is associated with 
the production of kinetic energy in the trailing vortex system, and is 
called "vortex drag." This drag is identical with that produced by the 
same spanwise lift distribution in an incompressible flow, (frequently 
called "induced drag"). 

The remainder of the drag due to lift and all of the drag due to 
thickness is associated with the production of energy near the surface 
of a downstream Mach cone whose vertex is in or near the aircraft. This 
Is called wave drag, and the associated energy is half kinetic and half 
( 1 ) 

potential' ' . 

The wave drag plus the vortex drag is equal to the drag evaluated 
at the wing and body surfaces by the first method, (it may be necessary 
to retain nonlinear terns in the expression for pressure coefficient to 
get this agreement . ) 

The momentum theorem is utilized in both of the above methods but 
different "control surfaces" are used. In the first case the control 
surface Is close to the aircraft surface, but in the second case the 

control surface is a distant one. For example Hayes uses a circular 
cylinder with axis passing through the aircraft and parallel to the free 
stream direction. The radius of the cylinder is chosen to be very large 
compared to the aircraft dimensions since this simplifies the ' calculations 

The wave drag of the aircraft is then evaluated from the rate at 
which momentum (in the free stream direction) is carried across the sur- 
face of the cylinder. (if the control surface had been chosen as a sur- 
face containing streamlines instead of a perfect cylinder, then the wave 
drag would have appeared as pressure on the streamline surface.) 

The cylindrical control surface is closed far downstream by a plane 
normal to the flow direction. The vortex drag is then determined, as in 
incompressible flow, by the rate at which the kinetic energy of the 
trailing vortex system passes through this plane, or alternatively 
through momentum and pressure considerations. 
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B^_ GENERAL MOMENTUM THEOREM FOR EVALUATION OF DRAG 

In the present section a momentum integral for the drag, as given 
by linearized theory, will be derived (Eqs. 4b-33,34). The drag will be 
given as an integral over an arbitrary control surface enclosing the 
solid. The integrand is a quadratic expression in the velocity compo- 
nents as given by linearized theory. 

First a more general momentum, integral will be considered. Consider 
a control surface S enclosing a solid (Fig. 4b-l). A surface element 
on S of area dS will be represented by its outward normal dn where^ the ^ 
length of dn is equal to the area of the surface element. Thus dn = (dS)n 
if ~n is the outward normal of unit length. Let the hydrodynamical stress 



Fig. 4b-l 


tensor be denoted by a, and let I be the region inside S and II the 
region outside S. Then 

f = a dn = Force exerted by II on X across surface element (4b-l) 

If a system of coordinates x^, X 2 , x^ is chosen dn may be repre- 
sented by its three components (dn) ^ and a by a 3 X 3 matrix ou j . The 
above equation may then be written 

3 

f ± = ( a dn) ± = ^ cr^dn^ (4b-2) 

J=L 

where (a dn)^ is the i™ 1 component of the force. 
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For a non-viscous fluid the only hydrodynamical force is the pres- 
sure p and the stress tensor is 

a = -pi = -(p8 it) ) (4b-3) 

where I is the identity tensor whose matrix is the Kronecker delta 5j y 
In this case the force across the element is 

f = -p(l dn) = -p dn (4b-4) 

or 

f ± = -p(dn) ± 


The hydrodynamical momentum equation states that the stress tensor 
is balanced by flow- of -momentum tensor. (This is actually a restatement 
of Newton* s law that force = (mass) times (acceleration).) To define 
the flow-of -momentum tensor we first introduce the concept of a dyadic 

product of two vectors. Let a and b be two vectors with components 
and The dyadic product is then the tensor whose ij component 

( a • b ) ^ j is b j j i . e • 


a . b = ^aj_ . b^ (4b-5) 

Note that if <c is any vector then 

(a . b )c* = ^z( a ibjjcj^ = ^aj_Ebjc^ = a(^ . c*) (kb-6) 

where b . c is the ordinary dot product. 

The flow of momentum tensor is then the dyadic product of pq 
(momentum per unit volume) and q velocity: 

^ ■> 
pq . q 


Flow-of -momentum tensor = 


(4b-7) 
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Its physical interpretation may be seen by applying this tensor to the 
normal dn 

(pq . q)dn = pq(q . dn) =VMomentum transported through dS per unit time . 

(4b- 8) 

The basic momentum equation for stationary flow for a surface S-j_ 
which does not enclose a body is then 




0 dn 


(4b-9) 


This is analogous to the law of conservation of mass which states that 



(4b-10) 


Consider now the composite surface consisting of the surface S in 

— 

Fig, 4b -1 and the body surface Z. Let dn denote normals on Z which point 
outwards with respect to the body (i.e. into region I). From the defini- 
tion of the stress tensor a 


F = Total force exerted by fluid on body 


■ X 


(4b-ll) 


Since the flow through S is zero one obtains by applying Eq. (4b-9) to 
the composite surface S^_ = S + £ 7 _ 


/ (pq • q)<in = / a dn - Fa dn 
J a J a J y 


(4b- 12) 


The minus sign in the last term is due to the convention that on the 

surface 2 the quantity dn denotes the inward normal with respect to the 
region I. Comparing Eqs. (4b-ll) and (4b-12) one obtains 
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F = - / (pq • q)dn + a dn (4b-13) 

J S J S 


This is the fundamental momentum formula which gives the total hydro - 
dynamical force on the solid as an integral over a control surface 
enclosing the solid . 

Note that in Eq. (4-b-ll), the force is given by an integral of the 
stresses on the body surface. This is the "close” point of view for 
evaluating the force. Eq. (4b-13) shows how the same force may be evalu- 
ated from the distant point of view. 

A slight modification of Eq. (4b-13) will be needed later. Denote 
the flow quantities at infinity as follows 

9; P> P> o' at infinity = U, p Q , p Q , a Q} respectively (4b-l4) 


The difference "between a flow quantity and its value at infinity will 
be denoted by a "prime . " Thus 


■», -> 
q. = q. 


- U, p' = P - P Q , p* = p - p D , a' = a - cr, 


(4b-15) 


From the continuity equation (Eq. 4b-10) it follows that 


-> 

U 



-> 

pq. 


dn = 



(pU . q)dn = 0 


(4b-l6a) 


F\irt hermore , since a Q = Constant 



dn = 0 


Subtracting Eqs. (4b-l6a, b) frcm Eq. (4b-13) one obtains 


-5» 

F 


- f (pq/ • %.)& + f 

UR J R 


a' 


dn 


(4b-l6b) 


(4b-lTa) 
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where, for a non -viscous fluid. 



dn 


(4b -17b) 


This is the fundamental momentum formula in terms of perturbation 

quantities. Note that the latter are not assumed to be small . 

The drag is the component of F in the free stream direction. We 
shall take this direction as the x-direction and use the following 
notation. 

U = U ± , q = (u,v,w), q' = ( u ' , V 1 , w 1 ) (4b-l8) 

where 


u = U + u'. 


- 1 _ 


= v, 


w' = w 


From Eq. (4b-17a) then follows the fundamental momentum formula for drag : 
Drag = F . i = - / pu'q . dn + i . j a' dn (4b-19) 

JR JR 


The momentum integrals may be further simplified for special choices 
of the control surface S, in particular by letting S recede to Infinity. 
However, we shall first derive an approximate form of the drag formula, 
valid within linearized theory. In the following section this linearized 
formula will then be specialized to a special Infinitely distant control 

surface (method of Hayes ( 1 )) . 


Inviscid Second-Order Drag - .i — 

It will be shown below that for a thin or slender body the largest 
contribution to the drag may be evaluated by an integral of a quadratic 
expression of the linearized perturbation velocities. It is usually 
stated that the drag is of second order. "However, it should be remembered 
that the values of the perturbation velocities are computed from first- 
order (linearized) theory. The result is - a formula for drag according 
to first-order theory. The term "second-order drag" refers to the fact 
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the integrand is quadratic in u 1 , v 1 and w 1 and hence of second order 
if u ! , v 1 and w 1 are themselves of first order. Furthermore, the second- 
order correction to u T , v 1 and w* will contribute nothing to the second- 
order expression for drag. The final formula is given by Eqs . (4b-33j34) 
and the reader interested only in the final result may skip the deriva- 
tion now presented below. 


We shall first assume non-viscous flow, so that the stress tensor 
is given by Eq. (4b-3). Furthermore, we shall assume that the solid 
is characterized by a parameter e, which is small, e.g. the fineness or 
thickness ratio. We shall furthermore assume that the flow quantities 
may be expressed by power series in e: 


u = U + eu-^ + + • • • 

v = ev*! + €^2 + • • • 

w s ew^ + • 

p = P o + ep x + £^2 + . . . 


(413-20) 


P = P Q + ep x + e P 2 + 


Such an expansion is valid at a distance from the body. It should 
be remembered, however, that in slender body theory, terms involving 
log e are of importance very near the body. 

The coefficients of e are the first order terms and are given by 

linearized theory. The coefficients of are the second order terms, 
etc. The lowest order term in the expression for the drag will now be 
found using the isentropic pres sure -density relation and Bernoulli ! s 
law. 


Fran isentropy it follows that density is a function of pressure 
alone . One defines 


(*), 


= a c 


constant entropy 


where a is the isentropic speed of sound. Then 


P = Po + -Vp - Po) + 


(4b-2l) 
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from which then follows 


p -Si. 

P 1 a 2 


Bernoulli's law may he written 

d(u^ + + w2) 


+ dP = 0 


where 


or 


- r s st 

" Jr. » 


(u 1 + U)^ + + w 2 + -2P = U c 


Using Eq. (4b-2l) P may he expanded to second order - 


*/■ 


dp ' 


P * 

P ° + a 0 2 


1_ 

P~ 


1 - 


dp' = — 


Po a o £ 


P' - 


(P ’ ) 


m2 


2p o a O^ 


Expanding the terms in Eq. (4b -25) to second order one obtains 


o P U l 2 + Y l 2 + V l 2 1 I o *=“^1“ \ 

eUu! + s%u 2 + c 2 -i 1 ^ + ^-(cPi + e2p 2 - = 0 


2pi 2 


2p 0 a o 


Collecting the terms of order e one finds the linearized Berno ull i ' 
law 


P 1 + Po u l U = 0 


(4b-22) 


(4h-23) 


(4b-24) 


(4b - 25 ) 
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Comparing with. Eg. (4b-22) one sees that 


Uu x 

Pl = " P ° 712 


The terms of order yield the following expression for p r 


(4b-26) 


Pq U u 2 + Po 


p P P 

U-^ + + Wj*- 


p Ui % 2 

+ ^2 - “ 2 4— = 0 


(te-27) 


where Eq. (4b -25) has been used and 


M = 


U _ 


In the momentum formula. Eg. (4b- 17 ), the stress and momentum flow 
tensors may he combined to form a tensor A 


A = - pg ’ . g - p'l 

Using Egs. (4b-20, 25, 26) one may evaluate An 


An = “P u ' (U + u' ) - p’ 


e (po u l U + Pl)" e2 (p 0 u 2 U + .Pl u l U + Po u l 2 + P2) 

/ u -, 2 + v , 2 + w -, 2 Ui % 2 N 

€ 2 Po (uu 2 - M 2 ^ 2 + U 1 2 - UU 2 - -i 1 i- + 


Finally 


An = e 2 P 0 § 


^M 2 - l)^ 2 + Vi 2 + w l2 ] 


(4b-28) 


Similarly 


A i2 = - pu ' v ’ A -€ 2 p 0 uivi 


(4b -29) 
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A 1J = -P u ' v ' = -e 2 P 0 u i w i . (413-30) 

Since only the first row ^A^, &12’ ^ 13 ) enters in the drag computation 
we have proved the following: 

1. The dominant term in the drag formula is ofl.second order in e 

2. The integrand in the drag formula is, to second order, a second 
degree polynomial in the first order velocity perturbations. The velocity 
perturbations of second order, or pressure and density perturbations of 
second order, do not enter into this expression. _ 

Thus while drag is of second order, it may be computed on the basis 
of first order theory (linearized theory) . On the other hand, one may 
easily check from the above expressions that in general the lift has a 
first order term. Furthermore to compute lift to second order one needs 
to know U 2 , that is, u to second order. 1 7J_ 

In the remainder of this report we shall only be concerned with the 
drag as given by linearized theory. It is then convenient to introduce 
a change of notation. We shall let u, v, w stand for the linearized 
velocity perturbation ; in other words 


eu^_, €V-]_, €W^_ are replaced by u, v, w 
Furthermore a velocity potential (j) will be introduced such that 


(4b-31) 


Grad <f = u, v, w 


(4b-32) 


The above results may then be summarized as follows. The drag to 
second order is given by the formula — - 


D = 



(4b-33) 


where S = Control surface enclosing the body 
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^1 = ( A ll> A 12> a 15 ) 

A 11 = + Po ^( p2u2 + v 2 + w2 )> H 2 = M 2 - 1 
A 12 = "Po uv 


(41-54) 


A 15 = -Po uw 

and u, v and w are the .components of the perturbation velocity given by 
linearized theory. 


C . HAYES METHOD FOR DRAG EVALUATION 


The method developed by Hayes in Ref. 1 consists in applying the 
drag formula Eq. (4b-35) to a special control surface, a truncated cir- 
cular cylinder, surrounding the body and in considering the limiting 
case when the control surface recedes to infinity. The general momentum 
integral for the drag then assumes a simplified form. (This results in 
certain simplifications in the integrand.) Furthermore, if the body is 
represented by singularities (sources, lifting elements, etc.) as dis- 
cussed in Ch. Ill, the velocities at large distances may be represented 
very simply in terms of the strength of the singularities. As a result 
the drag may also be represented as an integral over the singularities 
(distribution of source strength, etc.). This result of Hayes ' gener- 
alizes a previous result by von Karman(T) for a body of revolution. 

First a somewhat detailed demonstration of the method of Hayes will 
be given for the case of a lineal source distribution. This part may 
be skipped by a reader not interested in mathematical details. Then the 
results of Hayes and related results will be stated in intuitive terms 
for general three -dimensional distribution of sources, lifting elements 
and side-force elements. Detailed proofs will not be given. However, 
the results may be proved by methods closely analogous to the method 
exhibited for the case of a lineal source distribution. 


Hayes Method for Lineal Source Distribution 

We shall consider a distribution of sources along the x-axis between 
x = 0 and x = L. The corresponding solid is then a body of revolution. 
The source strength will be denoted by f . It will be assumed that 


f(0) = 0, f(L) = 0 


(4c-l) 



56 


NACA TO 1421 


These assumptions lead to certain restrictions on the body shape. 

Let the radius of the body be r(x). The cross sectional area S(x) is 

then «r2(x) . Since f(x) = U S'(x), f(0) means that r(o) . r'(o) = 0. 

This is fulfilled if r ~ x n , n > l/2 near the origin; - In particular, 
f (0) is equal to zero if the body starts in a point with finite slope, 
i.e. r ~ x near x = 0. The analogous condition at x = L insures f(L) = 0. 
In addition, f (L) = 0 if the body ends smoothly in a cylinder with con- 
stant radius, i.e. if S(x) = Constant for x > L and S'(x) is continuous 
and hence zero at x = L. It will be indicated in the proof below why 
the restrictions on f are necessary. 



Then t| = 1 on the downstream Mach cone from x = g and 0 < t < 1 inside 

this Mach cone. For x - pr ^ L one may also write the_ velocity components 
as 
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0x = f(i) [(x - S) 2 - P 2 ** 2 ^ (x - £)d£ 

= -^ f f (l)( x - O 2 (l - t| 2 ) d| (4c-5a) 

0 r = jT f(l) [(* - £ ) 2 - P^J p 2 r d| 

1 ~ [ 3 /2 

= jf f(0(x - - t| 2 ) d| (4c-5b) 

Hayes' Control Surface 

Following Hayes we now introduce the control surface shown in 
Fig. 4c -1. It consists of a circular cylinder of radius r^, truncated 

hy a front disc x = Constant < 0 and a rear disc x = x^_ > L. The drag 



Fig. 4c-l 
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integral (Eq. 4b-33) will be evaluated for this control surface as r^ 
and x^ tend to infinity. The ratio between x^_ and r-j_ will be determined 

later in such a way that the contribution of the rear disc to the drag 
will vanish in the limit. 


Contribution of Rear Disc _ 

According to Eqs. ( 1 rt)-33>3*0 the contribution of the rear disc to 
the drag is 

D = 1 (p^x 2 + ^r 2 ) 2 * 1 * dr (4c-6) 


The velocity components may be evaluated as follows. Write f(£) 
as a difference of two positive functions 


f(£) - f+(6) - f-(S), f+(6), f-(i) - 0 


_(4c-7) 


Then by the mean value theorem and Eqs. (4c-5a,b) 




f f+(ft)d£ 

V C\ 


f L f-(l)d£ 

C\ 


- 5 3 ) 2 (i - % f 2 *(* - rff - 


(4c-8) 


where 0 < £^, < L. A similar egression is valid for 0 r * Note that 

in Eq. (4e-8) the continuous source distribution is replaced by a posi- 
tive source at £j and a sink at £ 2 . However, £j and £ 2 depend on x 

and r. 

As is easily seen 


O^x-L^x - £^, 


t L ^t gj , i = 2, 3 


1 < 1 + tg , 1 + t io < 2 


'i- 

*2 


(4c— 9) 
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Hence, replacing and I 2 by L increases the absolute magnitude of 
both terms in Eq. (4c-8). Hence, on the rear disc x = x-^, 0 < v <, r^, 



where A is independent of Xj_ and r^_, and 



(4c -10) 


(4c-ll) 


If one puts y = 1 - tj^, then dy = dt^ = -t-^ . 2(3 dr^x-^ - 


Hence 


t 1 r3 2 ^ x 2 dr ^ 

on 


( X 1 " L ) 2 ^ yi y5 2(xi - L) ; 


Vfl 

y 2 Ji 


where y^ = 1 - ^1 - = 2e^_ - e^_ 2 and e-^ is explained in Fig. 4c-l. 

Equation (4c-12) may be written 


r 1 r|3 2 ?S x 2 clr < 


K-^1 2 r i 2 * 2 


(4o-13) 


where C is independent of r-^ and x-j_ for r^_ and sufficiently large; 
e is explained in Fig. 4c-l. The fact that *1, e / e l'~ > * i - 
as r-^>« has been used above. It is then clear that if e is constant 
or if e = r~ a , a < 1, then the integrand in Eq. (4c-15) tends to zero 


A similar estimate may be shown for /r0 r %x. A comparison 


with Eq. (4c-6) shows that: 
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Contribution of rear disc to drag is zero -even If e decreases 
as r increases. However, e should, decrease more slowly than r~^~. 
Since the distance BC is of the order er it follows that thiB 
distance becomes infinite in the limit. 


Contribution of Cylindrical Fart 1 . . _ 

Since the contribution of the forward disc to the drag integral is 
identically zero it follows then in the limit r-j_, x^_ — *■<» the entire drag 

contribution comes from the cylindrical part, provided e varies aB pre- 
scribed above. Thus 


D = Limit D 2 


(4c-l4) 


where Dg, the contribution of the cylindrical part, is 


ox+pr-^e 

d 2 “ -Po 2jac i / dx ' (4 c -15) 

J Pri - 

(Note that the radial component A-jj, of the vector Aj_ in the. drag formula 
Eq. (4b-34) is - p o 0 x 0 r • ) the above equation l/l - e has been replaced 
by 1 + e which may be done without loss of generality. 

To evaluate D 2 , we write Eq. (4c-Ja,b) in the following form 



0r = 
where x' 


-1 r x '> L 



0 

- e x f*' - ?! + 

2gr x 

p x ' } L 

f'(i a )|B 


2* J 0 /x 1 ” 

- t 2 i x ' - ^2 + 

2Pr 1 

* x - J3r-j_. 


_ 


+ - 
f*l ‘ 



(4c - 16 ) 
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The upper limit is x' for x* < L and L for x' ^ L. Hence 


nPeri ^x-l'L px-l'L x 

’ f '(Sl) f ’(52) 


o 

o 

o 

ro 

- ll l/x' - 12 

^ x 1 - | L + 2pr L 




x' + pr-L - £, 


fiT’ - | 2 + 2 P r l pr i 


d| 2 dx' 


(4c -IT) 


Limiting Case for Infinitely Distant Control Surface 

We shall now evaluate D 2 as — >». The three ratios within the 
second bracket all tend to unity as and may hence be neglected 
in the limit . Note that this approximation implies 


1 - 1 i- i .1 

^x - S -^ 2 - p 2 r 2 ^(x - ip) - pr |flpr ^x' - |fipr 


(4c-l8a) 


Furthe rmore , applying the same approximation to Eq. (4c-l6) one obtains 
that 




y x 


(4c - 18 b) 


0 r and 0 X both vanish as 1 j]j 2pr^ . Their ratio, however, is given by the 
above equation. The corresponding relation with 0 r replaced by is 

exact for two-dimensional flow. Thus the flow is approximately two- 
dimensional at large distances near the Mach cone from the leading edge 
small, i.e. t| almost unity for 0 ^ | ^ L^. 


Hence 


D = Limit D P = ?£. I 

* 4* 


(4c-19a) 


where 


operx r x',L r x’,L 
I = / • / / 7 aLa 7 === === dll dig dx- 

J 0 Jo Jo jfx* - I x |fx* - Ig 
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The domain of integration is- a region in x', £ 2 space whose cross- 

section for x 1 = Constant is the square 0 £ — X S 0 ^ £ 2 — x ' £° r 

0 ^ x' ^ L and the square 0 ^ l 2 — L f° r L — x’ ^ ’Per^. Let Lj_ 
be the integral where 0 < y < L and I 2 the integral over the domain 
L ^ y Per^. Since the integrand is symmetric in £]_ and | 2 , half its 
value is obtained by integrating only over the triangle ABC in the £^, 
i 2 plane (i.e. 1-^ — S 2 ) 613 s k own is. Fig. 4c-2. In evaluating over 



its domain (a truncated triangular cylinder with base at x = L and vertex 
at x' = 0) we shall first integrate along a line parallel to the y-axis. 
For £-±_, £ 2 fix®*! this line is inside the pyramid only when ^ x ' — 

may vary inside the triangle between 0 and | 2 , and for £2 an y value 
between 0 and L may be chosen. - Hence Ij_ may be written 


P L r i 2 p L 

I L = 2 / / / ^==S =- } ~ ks^—- dx' dE-L d| 2 

Jo ^0 J £ 2 P - ii - s 2 

The integral I 2 is (domain is triangular cylinder) 

1 . 2 r L r Ss d6 d , a*. 

J l J 0 Jo ^x' - ^ ^x’ - 1 2 1 2 


(4c-20) 


(4c-21) 
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Interchanging the order as above one obtains 


I -11+12.-2/ / f , (il)f , (l 2 )l / Cl d|i d| 2 

-J (4c-22) 

Here C is a constant and it has been introduced under the assumption that 


P f'(|)dg = 0, which, since f(0) =0 means f(L) = 0 
J n 


(cf. Eq. 4c-l) . ^Note that the limit of integration for may be 
replaced by L if the factor 2 is omitted.) 


Now 


L 


P er i 


.tel 


$2 - ll - S2 


- log C 


- |^s(2x’ - - 6 2 + 2 ^(x' - Ij) (x’ - i 2 ) 


per]_ 

12 


log C 


/ X 22^1 - +2lf/p€i* 1 -^V&er -fc 2 ) 

« -log(e 2 - ix) + log 1 g ' V c ■ ^ i— ^ (4 c-2J) 

f 

Hence if one chooses C = the second term will tend to log 1 = 0 

as r^ — Note that for this it is essential that er^ — * « as r-^ — 

(cf. pg.IVl6) . In other words the simplicity of the proof depends on 
the fact that € — >0 as rj _ — >*> (cf . Eqs. 4c-l8). On the other hand e 
may not tend to zero so fast that er^ remains hounded. In this case the 

above proof would he invalidated. Actually a drag contribution would 
come from the rear disc in that case. 

By combining the Eqs. l4-c-19a, 22.> 23 one obtains the final drag 
formula 


D 



^2 




(4c-24) 
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This is von Karman's drag formula for ~.~a lineal source distribution 
such, that f(0) = f(L) = 0. It has been derived above “by the method of 
Hayes . This derivation has the advantage that it may he extended 
Immediately to cases of a more general distribution of singularities . 
Such generalizations will now be discussed. 


General Three-Dimensional Source Distributions 


We shall now consider a more general case of a spatial distribution 
of sources. It will still be assumed that no lifting or side force ele- 
ments axe present. The source strength will be denoted by f(x,y, z). 

It will be assumed that f = 0 outside a certain finite region V. A 
special case is a planar distribution, say in the plane z = 0 in which 
case f(x,y,z) = f£(x,y)8(z). Another special case is the lineal distri- 
bution on the x-axis which was discussed above. In this case 
f(x,y,z) = f-j_(x)6(r). It will be shown below how in a certain sense the 

drag evaluation for the general three-dimensional case may be reduced 
to a consideration of certain equivalent lineal distributions. In the 
course of this discussion certain restrictions on f(x,y, z) will be made 
in addition to the requirement that it vanish outside a finite region. 

Consider a line in the streamwise direction passing through V. The 
position of the line which will be taken as the x-axis is actually arbi- 
trary, but for practical purposes it will be assumed that it is "well 
centered." This is, of course, a somewhat vague requirement. However, 
if for example f has rotational symmetry, the x-axis will be its axis 
of symmetry. On the x-axis choose as origin a point, 0, whose downstream 
Mach cone contains V. For convenience choose this point as far down- 
stream as possible. Also choose the point, L, for convenience as far 
upstream as possible, whose upstream Mach cone contains V. An equivalent 
requirement is that the downstream Mach cone from L is contained within 
the downstream Mach cones from every point in V. Let the value of x at 
point L be L. Thus the downstream Mach cone from x = 0 and the upstream 
Mach cone from x = L touch but do not penetrate V. We now introduce a 
control surface and define e and as in the lineal case (cf. Fig. 4c-l) 

This is shown in Fig. 4c-J. It will be assumed that r^_ and Xj_ tend to 
infinity as described in discussion of the lineal case. 
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Fig. 4c- 3: Hayes control surface in three-dimensional space 


That is as and r-^ tend to infinity e and will tend to zero. In 

that sense the line AC will come arbitrarily near the Mach cone from the 
origin. On the other hand e and €]_ will tend to zero slower than l/r^ 

so that the line AC becomes infinitely long as r^ — >oo. 

By the same methods that were used in the lineal case, it may be 
easily seen that the contribution of the rear disc, x = x-^, becomes zero 

in the limit. All the drag thus comes from a portion on the cylindrical 
surface arbitrarily near the Mach cone from the origin and is hence pure 
wave drag. 

To evaluate the drag contribution from the cylindrical surface we 
introduce cylindrical coordinates x, r, 0 where 

x = r cos 0, z = r sin 0 


(^c-25) 
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Let the drag contribution of a strip on the cylinder between 9 = 0 O 
and 0 = 0 O + A0 be AD. We define 


— = Drag contribution per unit angle = lim ^ as A0— >0 (4c-26a) 

d0 AS 


Then 


D = Total drag = 



(4c -26b) 


Consider now a fixed meridian plane 0 = 0 O , and a point P = (x^r-j^So) 

on the cylinder between A and C (Fig. 4c-5). The potential 0(P) depends 
on the contribution from all sources inside the upstream Mach cone from P. 
The contribution from a source at Q = (£ ,T], £) is proportional to the 
source strength f (Q) and inversely proportional to the hyperbolic dis- 
tance r^P^Q) between P and Q where - 

(i^cos 0 O - ti) 2 + ( ti _ sin 0 O - £) 2 1 (4 c-27) 


rh 2 = (*6 - 6) 2 - P 2 


This hyperbolic distance is constant on hyperboloids of. revolution with 
r = 9 = 0 O as axis . Consider now the sources between two such hyper- 

boloids which intersect the x-axis at x = | and x = £ + d£. To evaluate 
the contribution to 0(P) of these sources one may transfer their total 
source strength to the axis. In this way the distribution in V is 
replaced by an equivalent lineal distribution i.e. by an equivalent body 
of revolution. So far this lineal distribution depends on x Q and r^ as 

well as 0 O . 


Consider now, still for fixed 0‘ = 0 O , the limit as r-^ — Then 

the hyperboloids may be replaced by Mach planes which intersect the 
meridian plane 0 = 0 O orthogonally along Mach lines. .Note that for this 

it is necessary that as — >°° any point between A andj C comes arbi- 

trarily near the downstream Mach cone from the origin .in the sense 
described above. The source strength between two such neighboring planes 
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Fig. 4-C-4-: Evaluation of 0(P) 

may then he transferred to the x-axis as above. However, in this limiting 
case the resulting equivalent body of revolution depends at most on 0 O . 

It becomes independent of r^ and Xq. The corresponding lineal source 
distribution will be denoted by f (x$9o) • A consequence of the independ- 
ence of x Q and r-j_ is that f (x;9 ) may be used for computing 0 r and 0 X as 
well as $ at P. In general it may not be 'used for computing 09 . Clearly 
0 q is zero for a lineal distribution, whereas the 00 resulting from the 
original volume distribution is not. On the other hand 0 0 is not needed 
for drag evaluation on the cylindrical surface. 

Since 0 r and 0 X may be computed from the equivalent body of revolu- 
tion for fixed 9 it follows that <3D/d9 may be computed in exactly the 
same way as the drag of a body of revolution was computed. The result 
will differ from Eq.. (4-C-24) only by a factor 2it. Hence we have proved 
the following: The drag D of a volume distribution of sources of 

strength f(|,rj,^) is given by the formulas 
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D = 



££ d.9 
de 


(4c - 28 a) 


dD 

d0 



(^2^ ®) ^-°S (^2 “ ^2 


(4c - 28 b) 


f(t;Q)ai 


-Iff 

V(ft,0) 


f (Q)dQ 


(4c - 28 c) 


where V(g,0 o ) is the region contained "between two Mach planes perpen- 
dicular to 0 = 0 O and intersecting the x-axis at x = g and x = g + dg . 

This result was obtained by Hayes in Ref. 1. It is thus seen how 
Hayes' derivation of von Karmen's drag formula for bodies. of revolution 
admits an easy generalization to the general three-dimensional case. 

This proof obviously presupposes the following' requirement on the 
strength distribution f(Q) in addition to the requirement that it vanish 
outside a finite volume: f (Q,) must be such that for each 0 f (x;9 ) sat- 

isfies the same requirements as f(x) in the lineal case. In particular 
for each 0: f(O;0) = f(L;0) = 0 and f(x;0) must be differentiable with 

respect to x. 

If f (Q) has rotational symmetry, i.e. depends on r and x only then 
it may obviously be replaced by one equivalent lineal distribution, 
independent of 0 , for computing the distant flow field and the drag. 

In the special case when f(Q) is lineal to begin with, Eqs. (4c-28) 
reduce to the previously established formula (Eq. 4c-24). 


Extension to Include Lift and Side Force Elements 


For simplicity only sources have been considered in the preceding 
development. However lift and side force elements can be included and 
were included by Hayes in his original report. We will not go into the 
details here, but merely indicate the final results, since the funda- 
mental ideas of the method have been illustrated in the discussion of 
source distributions. 

Following Hayes we define a function h such that 

h = f - ^g z sin 0 + gy cos 0) (4c-29) 
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where 


f = f ( | ; 0 ) = Source strength 


pUg z = l(|;0) = Lifting element strength 


pUg = s(|;0) = Side force strength 


The term ^gy sin 0 + g x cos QJ is proportional to the component of force 

in the direction 0, and is the only component contributing to the wave 
drag in the Hayes calculation. Equation (4c-28b), as extended to include 
lift and side force elements , is 


Slave - Sri 5 ' - 6 l) dS l «i 


-p ,-L P L 
J 0 J 0 


f f h'(i i; e)h'(i 2 ;e)io E 


h - d s 2 


(4c - 50 ) 


where h(|;0) is the equivalent lineal distribution (for a given station 0) 
of the original spatial distribution of singularities. 

This equation makes it possible to determine the wave drag of an 
arbitrary spatial system containing thickness and carrying both lift 
and side forces. In order to determine the total pressure drag of the 
system it is necessary to evaluate the vortex drag produced by the lift 
and side force . In Hayes method the vortex drag appears as a mcmentum 
outflow through and a pressure on the end of the cylindrical control 
surface. It can be evaluated by calculating this momentum and pressure 
or by determining the kinetic energy associated with the vortex system 
in the Trefftz plane. Since this is identical with the induced drag 
problem of incompressible flow, we will not discuss it further. 


D. LEADING EDGE SUCTION 


The evaluation of the drag of a lifting wing of zero thickness by 
integrating local pressure times frontal area over the wing surface is 
not .theoretically complete until leading edge suction is accounted for. 
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This means that the infinite negative pressures acting on subsonic 
leading edges should be included. In practical applications this leading 
edge suction is sometimes discarded since in many cases only a fraction 
of the theoretical value is actually realized. 

However, from the distant viewpoint, leading edge suction cannot 
be isolated. This is true because there is no point-to-point corre- 
spondence between the close and the distant control surfaces. At the. 
distant control surface the velocity field created by the wing leading 
edges is merged with the fields created by other areas on the wing and 
body . 


From the distant point of view leading edge suction is automatically 
assumed to be fully effective, and therefore it must be so assumed from 
the close viewpoint to get correspondence in the drag values. 


E. DISCONTINUITIES IN LOADINGS 


For a planar wing, vortex drag is dependent only on the spanwise 
lift distribution. A discontinuity in the ordinates of this lift dis- 
tribution produces a concentrated vortex of finite strength and infinite 
energy, which corresponds to infinite drag. 

Wave drag is similarly affected by discontinuities in loadings. 

For example, consider a distribution of sources on a streamwise line. 

If there Is a discontinuity in source strength, then the drag evaluated 
on the distant control surface is Infinite . 

To prove this, assume a source distribution with a discontinuity 
at the point x = x (see sketch). The velocity potential at a point (x,r) 
downstream of the rearward Mach cone from x may be written 


0 


■ ' 


x-j3r 




JC 


x-pr 


1 /(x - \ ) 2 - p 2 r 2 

£g(Od| 

^(x - £) 2 - p 2 r 2 


JL_ 

2jc 


L 


X 




o ^(x - 1)2 _ p2r2 


(4e-l) 


-f(f) 



o 


X 
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The u-component of velocity at the point (x,r ) is found hy differentiating 
Eq. (4e-l) axi ally , (in order to avoid indeterminant forms in the differ- 
entiation, the equation is first transformed by means of the relation 
6 = x - pr cosh u.) This process gives the result (assuming f(0) = 0): 


M 

n x f^g) d| nx-pr f 2 «(|)d| 

Af(x) 

dx 2rt 

~'0 |/(x - | ) 2 - p^ 2 Jx |/(x - |) 2 - p^ 2 

^(x -x) 2 - p 2 r 2 


where Af(x) = f]_(x) - f 2 (x) . 


(Jte-2) 


At the distant control surface it previously was shown (Ch. TV-C) 
that one need consider only conditions very near the Mach cones from the 
source distribution. Introducing the approximations used in Hayes' 
method (i.e., (x - |)/pr 3* l), Eq.. (4e-2) can be expressed 


u 



Kx' - l 


+ 



' f 2 ,(S)d| 
U' - l 


Af(x) t 
Px’ - x 


(4e-3) 


where x 1 = x - pr and x‘ - | « pr. Since the radius of the control sur- 
face is large compared to the length of the source distribution, the Mach 
cones originating at the sources are essentially plane waves when they 
intersect the control surface, so that the radial component of velocity 
(at the control surface) is (Eq. 4c-l8b) 


v = pu 


(4e-4) 


The drag, being equal to the transport of horizontal momentum across 
the control surface, is proportional to the product of u and v integrated 
axially along the control surface. Fran. Eqs. (i+e-3) and (4e-l+) it is 
readily seen that the drag includes a term of the form 



dx' 


The integral is non- convergent . An infinite drag contribution therefore 
results from a discontinuity in the strength of the source distribution. 
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F. THE USE OF SLENDER BODY THEORY WITH THE DISTANT VIEWPOINT 


If slender body theory Is applied, then the source strength Is 
assumed proportional to the rate of change of cross-sectional area, 
dS/dx, for a corresponding body of revolution. This means that Infinite 
drag will be predicted (by the distant procedure) for all bodies of 
revolution having discontinuities in d£>/dx. Such a prediction is, of 
course, incorrect, and the error is caused by the application of slender 
body theory to bodies which are not sufficiently smooth. 

The use of slender body theory requires that smoothness should be 
maintained at the nose and tail of the body and therefore dS/dx should 
be zero at these locations. In order that dS/dx should be zero at the 
nose or tail of a closed body of revolution it is necessary that the 
variation of body radius, R, with distance, d, from the nose or tail 

should be of the form R ~ where k > 0. This does not elimi- 

nate blunt noses or tails entirely, but excludes "excessive" bluntness, 
(Note that the Sears -Haack optimum shape is blunt.) 

The linearized theory requirement that all velocity perturbations . . 
be small theoretically excludes al 1 bluntness, but this is unimportant 
if very small regions of the flow field are affected. 

Bodies which begin or end in cylinders also may satisfy the smooth- 
ness requirements. “ 

For a body to be sufficiently smooth to permit the use of slender 
body theory, it is necessary to restrict the "short" wave length fluctu- 
ations in the plot - of cross-sectional area versus length. ( The word 
"short" cannot be defined exactly here, but should probably' apply to 

all wave lengths less than the body diameter times |/m^ -~1J 

Figure 4f-l illustrates the effect of wave length on the accuracy 
of the slender body theory. The drag for an infinitely long corrugated 

cylinder according to strict linear theory was found by von Karman^). 
Slender body theory is in good agreement with these results only where 
the reduced wave lengths are large compared _to the cylinder radius. At 
the other extreme two-dimensional theory is approached. 

It should be remembered that when the distant viewpoint is used 
the drag of a singularity distribution is evaluated. The body shape 
corresponding to the singularities may be determined either by "exact" 
linear theory or approximated by slender body theory. For example in 
Fig. 4f-2 a specific source distribution is considered, and is inter- 
preted as a "bump" on a cylinder by "exact" linear theory and by the 
slender body approximation. For this ratio of wave length to cylinder 
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COMPARISON OF THEORETICAL CALCULATIONS 
FOR DRAG OF CORRUGATED CYLINDER 


(A) RATIO OF THE DRAG COMPUTED BY SLENDER BODY THEORY 
TO THE DRAG COMPUTED BY L /HE A R THEORY 
(B.) RATIO OR THE DRAG COMPUTED BY TWO -D/M ENStONA L 
UNEAR THEORY (ASSUMING TWO -DIMENS/ONA L FLOW /N 
EACH MERID/AN PLANE) TO THE DRAG COMPUTED 
BY THREE-DIMENSIONAL LINEAR THEORY 



(FOR L/NEAR THEORY DRAG OF CORRUGATED CYLINDER, SEE REF: 7) 


Fig. 4f-l 




vn 
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COMPARISON OF BUMP SHAPES OBTAINED FROM A GIVEN SOURCE DISTRIBUTION 
ACCORDING TO SLENDER BODY THEORY AND UNEARiZED THEORY 


U. 


SOURCE DISTRf BUT/ON /(x) 




f(x) - -o.tz TrLfoZsli-^f 

o ’o 

(rws SOURCE DfSTR/BUT/OK COULO BE USED 
TO FORM (ALONG THE AXIS') A SEARS- HA A CK 
BODY OF FINENESS RATIO S, LENGTH JOz) 
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diameter tbs "bump shapes and locations are quite different. It is of 
interest however that the net volumes contained in the bumps are iden- 
tical. This has been proved by Lagerstrom and Bleviss and generalized 
by Bleviss in Ref. 22. (This suggests that "volume elements" may retain 
their significance even when slender body theory does not apply. ) 


CL THE DEPENDENCE OF DRAG- COEFFICIENT ON MACH NUMBER 

Hayes^) has pointed out that, for a distribution of singularities 
on a single streamwise line, the drag, evaluated from the distant view- 
point, is independent of Mach number. If the singularities are sources, 
and slender body theory is applied, this indicates that the drag of a 
given body of revolution is independent of Mach number. However the 
application of slender body theory in conjunction with the distant view- 
point requires that dS/dx = 0 at the tail of the body. 

Hayes' result is therefore consistent with a fact previously deter- 
mined, that the drag coefficient of a slender body satisfying the 
"closure" condition (dS/dx = 0 at the tail) is independent of Mach 
number . 

If the singularities are not confined to a single streamwise line, 
then the distant viewpoint gives a drag coefficient which varies with 
Mach number. This can be seen from the fact that the projection of the 
singularity distribution onto a single streamwise line varies with the 
inclination of the Mach planes used for the projection. 


H. SUPERPOSITION PROCEDURES AND INTERFERENCE DRAG 


In all the developments discussed in this report the linearized 
supersonic flow equation is vised. This means that one flow field and 
the lift (or volume) distribution which causes it can be 'superimposed 
on a second flow field with its corresponding lift (or volume) distri- 
bution. If the individual flow fields satisfy the linearized flow 
equation, then their sum does also. 


For example, let a pressure field, p^, correspond to a downwash 
field, o^, and a second pressure field, p 2 , correspond to a second down- 
wash field, ct 2 , then the pressure field p-j_ + p 2 corresponds to the down- 
wash field a^_ + 02 - 


However, the drag of the sum of the two fields is not in general 
the sum of the drags of the individual fields. For example, the drag 


of the first field would be Dj_ = 



dS, 


where the integration extends 
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over the wing and body surfaces, and similarly the drag of the second 
field is I >2 = jT dB . However, the drag of the combination is 

D-j _ + 2 = Pi + P2)( a l + The terms involving cross products give 

the interference drag, + P^l)^- 


1^ ORTHOGONAL DISTRIBUTIONS AND DRAG REACTION PROCEDURES 

If the interference drag is zero then the two distributions are 
said to be orthogonal. The use of orthogonal distributions for reducing 
drag has been studied in Refs. 8, 9, 10, and 11. 

For example consider two types of lift distributions which are 

orthogonal and assume that each one carrie_s a net lift. It has been,, 
shown (see for example Ref. 9) that some combination' of the two will 
carry a given total lift with less drag than would be produced if either 
one of the individual types of distribution carried all of the lift. 

On the other hand, any given (non-optimum) lift_distribution can 
be improved by adding the proper amount of a non- orthogonal type of dis- 
tribution which carries zero net lift. The improvement is obtained by 
utilizing negative interference drag. This can be seen as follows. The 
total drag of the combination is the sum of the individual drags plus 
the interference drag. The interference drag can always be made negative 
by proper choice of the sign of the distribution that carries zero nejt 
lift.. Also, since the strength of the zero lift distribution enters 
linearly into the interference drag, but enters quadratically into its 
individual drag, the magnitude can be so chosen that" the interference 
drag dominates. Thus the total drag of the combination can be made less 
than the drag of the given (non-optimum) lift distribution. 


J. THE PHYSICAL SIGNIFICANCE OF INTERFERENCE DRAG 


It has been stated that the interference drag, D^, 


ib 


/( Pl a 2 + P 2 a i^ 


dS 


where the subscripts designate the two flow fields which have been super- 
imposed, and the integration is to be carried over all surfaces . Assume 
that both flow fields are produced by thickness distributions. Then the 
a values are the body surface inclinations"”which correspond to dS/dx, 

the rate of change of cross-sectional area for the body. The /p -jOg dS 
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gives the drag produced by the pressure field of the first body acting 
on the cross-sectional area distribution of the second. The term ^P 2 «i dS 
has a similar interpretation . 

Assume that both flow fields are produced by lift distributions. 

Then J' p 1 a 2 dS is the drag created by the downwash field of the second 

distribution acting on the lifting elements of the first distribution. 

(The surface which supports the lift corresponding to p 1 must be inclined 

further because of the downwash due to p 2 *) 

Let the first field be produced by a lift distribution and the second 

by a thickness distribution (a body). Then J~ p^cc^ dS is the drag produced 

by the downwash field of the thickness distribution acting on the lift 
elements plus the drag caused by the pressure field of the lift distri- 
bution acting on the cross-sectional area distribution of the body. The 




P 20-1 dS gives no contribution to the drag in thiB case . 

Assume that the first field is produced by a lift distribution and 
the second by a side force distribution. The J' pja^ dS is drag corre- 
sponding to the downwash field of the side force distribution acting on 
the lift elements , while the J' P 2 a ]_ dS is produced by the side wash field 
of the lift elements acting on the side force distribution. 


K. II^rogERMCE AMONG LIFT, THICKNESS, AND SIDE FORCE DISTRIBUTIONS 

For planar distributions of lift and thickness (the lift being normal 
to the plane) there are no interference drag terms, and the two problems 
can be studied independently. However, for spatial distributions, inter- 
ference generally exists. This has been discussed by Hayes, and the 
physical meaning of the interference drag has been discussed in the 
preceding sections . 
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Suppose that a source and a 
lifting element are located as shown 
in Fig. 1+k-l, the direction of flow 
"being perpendicular to the page. 

Then the component of the lift which 
lies in the line connecting the two 
singularities causes all of the 
interference . If the lift element 
were located on the y-axis (corre- 
sponding to a planar wing problem) 
there would be no interference. 

For lift and side force ele- 
ment s, as shown in Fig. 4k- 2, there 
is interference between the force 
components which lie in the line 
connecting the singularities, and 
also interference between the com- 
ponents normal to the connecting 
line . 

If the side force element lies 
either on the y-axis or on the z-axis 
(as shown in Fig. 4k- Ja and b), then 
there is no interference. This can 
also be seen from symmetry considera- 
tions, which show that the lift ele- 
ment produces no sidewash at the 
side force element and similarly 
the side force element produces no 
downwash at the lift element. 


“ . 2 - 



Fig. 4k-Ja 

2 


if 


Fig. 4k- 5b 
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L. REDUCTION OF DRAG DUE TO LIFT BY ADDITION OF A THICKNESS DISTRIBUTION 

Consider tbs two-dimensional system sketcbsd in Fig. 4l-l. The 
cross-hatched area is a thickness distribution lying partly in the pres- 
sure field of a flat-plate wing. The relative geometry of the thickness 
distribution and the lifting surface are indicated in the figure . Also, 
the pressure distributions, relative to the two-dimensional pressure 2aq/p, 
sure shown in parentheses. 



As long as the pressure field of the thickness distribution does 
not intersect the flat-plate, the lift of the system is the same as for 
the flat-plate by itself. On the other hand, the interference between 
the pressure field of the flat-plate and the thickness distribution pro- 
duces a negative drag contribution, so that the total drag of the system 
(omitting friction) is 12-l/2 percent less than the drag of the flat- 
plate alone. Thus, the total lift in this case Is unaffected by intro- 
duction of the thickness distribution and a drag reduction is obtained. 

This example is related to the Busemann biplane. The result obtained 
illustrate s the fact that, in the general case (non -planar systems), 
sources and lifting elements have an interference drag. 
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CHAPTER V. THE CRITERIA FOR DETERMINING OPTIMUM DISTRIBUTIONS 
OF LIFT OR VOLUME ELEMENTS ALONE . 


THE "COMBINED FLOW FIELD" CONCEPT 

The idea of the "combined, flow field.", was introduced, by Munk^^ 

and extended by R. T. Jones^ 1 ^* 1 ^ . Consider a distribution of lifting 
elements in a free stream of given velocity. A certain downwash velocity 
and pressure are produced at each point in the field. If the direction 
of the free stream is now reversed without moving the lift elements" or 
altering the directions and magnitude of these lift contributions then 
in general different downwash velocities and pressures axe produced at 
each point in the field. ^ 

One -half the sum of the downwash velocities produced at a given 
point in the forward and reverse flows is called the jlownwash velocity 
of the combined flow field at that point. — A similar _ definition applies 
to sidewash velocity. One-half the difference of the pressures in the 
forward and reverse flows is called the pressure in the combined flow 
field. These definitions follow from the super-position of the 
perturbation velocity fields for forward and reverse 'flow. It should 
be remembered that in the flow reversal the lift distribution ( not the 
wing geometry) is fixed. 

The same ideas may be applied if other singularities such as sources, 
side force elements and volume elements are considered. When sources are 
used the signs must be reversed when the flow direction is reversed. A 
source in forward flow becomes a sink in reverse flow. 


COMBINED FLOW FIELD.. CRITERION FOR IDENTIFYING 
OPTIMUM LIFT DISTRIBUTIONS 


A necessary and sufficient condition for minimum wave plus vortex 

drag was given by R. T. Jones^^) in connection with^planar systems. The 
condition is that the downwash in the combined flow field shall be con- “ 
stant at all points of the planform. This result depends on the fact 
that a pair of lifting elements ha6 the same drag in" forward, and reverse 
flow, which is also true when the lifting .elements are not in the' same 
horizontal plane . Hence the above criterion can be extended immediately 
to lift distributions in space by requiring constant downwash (in the 
combined flow field) throughout the space. 
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TEE COMBINED FLOW FIELD CRITERION FOR 
IDENTIFYING OPTIMUM VOLUME DISTRIBUTIONS 


A necessary and sufficient condition for minimum wave drag due to 

thickness was given by R. T. Jones^ 1 *^ in connection with planar systems. 
If total volume is fixed then the optimum distribution of volume gives 
a pressure gradient in the combined flow field which is constant over 
the planform. 

As in the case of lifting elements this criterion can be extended 
to cover thickness distributions in space. It is then necessary for the 
pressure gradient in the combined flow field to be constant throughout 
the space. 


D. UNIFORM DOWNWASE CRITERION FOR MINIMUM VORTEX DRAG 


A necessary and sufficient condition for vortex drag alone to be 
a minimum is that the downwash velocity throughout the wake of the wing 
system shall be constant in the Trefftz plane. (The wake cross-section 
Is the projection of the wing system on the Trefftz plane.) This condi - 

tion was given by Munk v . 

If the wake of the wing system has an elliptical cross-section then 
a constant intensity of lift over the cross-section satisfies the above 
condition and gives the minimum possible vortex drag. (See Appendix V-l). 
In particular when the cross-section of the wing wake degenerates into 
a horizontal line, (corresponding to a planar wing) the familiar require- 
ment of elliptic spanwise load, distribution is obtained. 


E. ELLIPTICAL LOADING CRITERION FOR MINIMUM WAVE DRAG DUE TO LIFT 


In special cases elliptic loadings identify minimum drag configura- 
tions, as has been shown by Jones . Let the space containing the 
lifting elements be cut by a series of parallel planes each inclined at 
the Mach angle to the flow axis. Consider all the lift intensity cut by 
any one plane to be located at the intersection of the plane with the 
flow axis. If the resulting load distribution on the axis is elliptical, 
and if this is true for all possible sets of parallel planes (inclined 
at the Mach angle), then the wave drag is a minimum. 

In Eayes^) procedure for calculating drag (see Ch. IV) this con- 
dition corresponds to obtaining the minimum possible drag contribution 
at every angular position on the cylindrical control surface. 
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Such minima cannot "be attained in general since the condition is 
sufficient hut not necessary. However if they are attained and if the 
vortex drag is also a minimum then the more general criterion (constant 
downwash in the combined flow field) is satisfied. 


Fj_ THE "ELLIPTICAL LOADING CUBED" CRITERION FOR 

MINIMUM WAVE DRAG DUE TO A FIXED TOTAL VOLUME ~ 

Sears and Eaack^ 1 ^) in determining optimum shapes for bodies 
of revolution in supersonic flow have also determined sufficient condi- 
tions for identifying optimum distributions of volume elements within 
a prescribed space. 

We consider a distribution of volume elements within a prescribed 
space and ask how these elements should be arranged in order that they 
should cause the least wave drag while providing a fixed total volume. 

If the equivalent body of revolution for a given angular position 0^ 

on the distant control surface (see Ch. IV) conforms to the Sears-Haack 
optimum shape then the wave drag contribution at 9^ is a minimum. There- 
fore if the equivalent bodies of revolution for al 1 values of 0 are 
optimum shapes the total wave drag is a minimum. 

The density of the lineal distribution of volume elements repre- 
senting the Sears-Haack optimum shape corresponds to the cube of an 
elliptical distribution over the length of_the line.”. Hence if all the 
equivalent lineal distributions have this form an optimum is ensured. 

Such minima cannot be attained in general since the "Elliptical 
Loading Cubed" criterion is a sufficient, but not a necessary condition 
for minimum drag. When such minima are attained the more general cri- 
terion (constant pressure gradient in the combined flow field) is also 
satisfied. 


CL COMPATABILITY OF. MINIMUM WAVE PLUS VORTEX DRAG 
WITH MINIMUM WAVE OR MINIMUM VORTEX DRAG 


It is possible for minimum wave plus vortex drag to be obtained 
when neither the wave nor the vortex drag is individually a minimum. 

For example consider that the "space" within which lifting elements 
may be distributed is the planform shown in the figure. For the vortex 
drag to be a minimum it is necessary to maintain an elliptic spanwise- 
loading over b. This requires a finite load on "a" which in turn pro- 
duces infinite wave drag if the chord for "a" goes to zero. However 
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the minimum drag due to lift for the planform is certainly finite (load 
the end pieces only and consider them as isolated wings) hence minimum 
vortex drag is not consistent with minimum total, drag in this case. 

On the other hand, for a planar wing of elliptical planform minimum 
wave drag and minimum vortex drag are obtained with the same (constant 
intensity) lift distribution. 


H. ORTHOGONAL LOADING CRITERIA 


Optimum distributions can be identified also through orthogonality 

considerations 9) . The optimum distribution of lifting elements in a 
space is orthogonal, to every distribution carrying zero net lift and is 
not orthogonal to any other distributions. 

A similar statement can be made for the optimum distribution of 
volume elements alone (assuming for the moment that negative local vol- 
umes are not excluded). However if lifting (and side force) elements 
are introduced in addition to volume elements, then the criterion must 
be modified. For example the rotationally symmetric wing plus central 
body having zero wave drag is orthogonal to all singularity distributions 
although it contains a net volume.* 

The criteria discussed in preceding sections of this chapter have 
not been thoroughly investigated for cases involving lift and volume 
elements simultaneously. However, some material on interference between 
lift and volume distributions is given in Ch. IX. 


See p. 105 ff . Since the wave drag is zero the disturbances on a 
distant control cylinder are identically zero. Hence its interference 
with any other singularity distribution is zero. 
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APPENDIX V 


DISTRIBUTION OF LIFT IN A TRANSVERSE 
PLANE FOR MINIMUM VORTEX DRAG 

As stated by Munk's Stagger Theorem^), the vortex drag of a spa- 
tial wing system is not changed if all lift and side force elements in 
the system are projected onto a single plane normal to the flight direc- 
tion (see Fig. A5-1). Furthermore, if there are no side force elements. 



Fig. A5-1 - 


then Munk's criterion for minimum vortex drag is that in the Trefftz 
plane, the downwash in the wake must be constant. (The wake cross- 
section is defined as the projection of the wing system on the Trefftz 
plane.) Assume that the downwash' field associated with the optimum lift 
distribution is w = -w D and that a uniform field w = +w Q is superimposed 
on the original field in the Trefftz plane f then the resulting two- 
dimensional flow pattern is equivalent to a uniform flow around a solid 
body. Munk gives the expression for the lift distribution in the trans- 
verse plane in terms of the velocity potential of this new flow for 
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certain todies symmetrical with respect to the x-z plane; for example, 
if <f is the two-dimensional potential flow around an elliptic cylinder, 


then 



boundary 


and 


dvortex^in - 




Fig. A5-2 


where l and d are the lift and drag intensities per unit area in the 
transverse plane. For an ellipse oriented as in Fig. A5-2, the 

potential is^ 1 ®^ 


0 = w Q (a + b)cosh (k - ! 0 )sin i\ 


where 


y + iz 


fa 2 - b 2 cosh(| + iT]) 


The curve 6 = | Q corresponds to the boundary of the lift distribution 
in the transverse plane. From the above equations one obtains 


Z = 2 ot/M = 2pUw 0 (a + b) 

dz /e=g 13 

so that the lift intensity in the transverse plane must be constant to 
obtain minimum vortex drag. With S = «ab, the drag is 

( v o\r _ L 2 

Dvort eXhlin = ^lfj L " 4qS(l + a/b) 

where L is the total lift generated. Thus to obtain minimum vortex drag 
for a spatial distribution of lift whose Trefftz plane projection is an 
ellipse with one axis vertical, the lift should be distributed so as to 
give a constant intensity when projected on the Trefftz plane . 
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This proof can he extended to cases in which the projected lift 
distribution covers a rolled ellipse, as shown in Fig. A5-3- If only 

lift (and no sidef'orce) elements are 
allowed, Munk's criterion of con- 



stant downwash still holds, but the 
lack of symmetry precludes use of the 
formulas given above. However, the 
optimum lift distribution can be 
determined by a superposition of two 
symmetrical optimum distributions, 
as shown in Fig. A5-4. L-j_ and Lg 



Fig. A5-3 



Co) (b) (c) 


Fig. A5-4 


are constant intensity lift distributions over the elliptic areas which 
produce constant downwashes w-j_ and Wg over those areas. Because the 

governing equation is the Laplace equation, which is linear, the lift 
distributions Lq and L 2 and the flow fields they produce can be super- 
imposed. If = L cos 0 and Lg = L sin 0 and Fig. A5-4c is rotated 

through the angle 0, then Fig. A5-4c corresponds to Fig. A5-3* There 
is a uniform downwash w corresponding to the uniform lift L. Thus 
Munk's criterion is satisfied and the drag is a minimum. It can be 
shown by symmetry that the total interference drag between the lift 
distributions Iq_ and Lg is zero so that the drag of L is obtained 

simply by adding the drags of L^ and L g ; that is 
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kj 2 •*~'2 2 _ L^(a sla.^0 + b cos^#) 

D ' rorte %in ~ 4 q s(l + a/b) + 4qS(l + b/a) ~ 4q.S(a + b) 


It should "be noted that for this optimum rolled ellipse case there is 
also a uniform sidewash generated. If a distribution of side force 
elements were available, it would be possible to utilize the uniform 
sidewash to reduce the vortex drag below the value given above. 
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CHAPTER VI. THE OPTIMUM DISTRIBUTION OF LIFTING ELEMENTS ALONE 


THE OPTIMUM DISTRIBUTION OF LIFT THROUGH A SPERICAL SPACE 

Consider a sphere of radius "R" with its center at the origin, and 
let a total lift "L" he distributed through the sphere with local inten- 
sity "2." If 2 = - , r being the radial distance from the 

*%2/r 2 _ r 2 

origin, then elliptic loadings are obtained when the sphere is cut by 
any set of parallel planes (see Appendix VI for derivation). The fact 
that elliptic loadings are produced when the planes are inclined at the 
Mach angle (to the free stream direction) insures that the wave drag Is 
a minimum (Ch. V). The cross-section of the wake is circular, and. If 
the lift intensity is projected onto a plane normal to the free stream 
direction it can be shown that the lift is uniformly distributed over 
this circular cross section. This insures that the vortex drag is also 
a minimum (Ch. V). 

The lift distribution 2 = y ■ then gives the minimum pos- 

jt 2 R 2 ^R 2 _ r 2 

sible wave and vortex drag. By Hayes' procedure it can be found that 

l 2 b 2 

the minimum vortex 

2*q(2R)% 2 


and the minimum total drag Is 


The largest planar wing of circular planform contained in the sphere 

has a minimum drag^ 1 * 4 ') which i6 greater by the ratio — — . This is 

2M 2 - 1 

a factor of 1.885 at M = { 2 . However, the drag comparison is, of course, 
not complete without consideration of the viscous drag (and thickness 
drag). For the spatial lift distribution described above, the required 
wing area is infinite and so, then, is the viscous drag. But the same 
minimum of wave and vortex drag can be achieved with a number of wing 
systems having finite wing area. For example, consider_the Infinite set 
of cascades enclosed in a spherical space as shown in Fig. 6a-l. At 


the minimum wave drag is -wave 


drag ^5) is r^ in vortex = L 


2nq(2R) £ 


J min ~ 


1 L* 
*q(2R) ; 


2M 2 - 1 
M 2 
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Fig. 6a-l: Cross-sectional view of an optimum set of finite area 

lifting surfaces in a spherical space 


M = {2 this set of cascades covers the region adequately so that the 
equivalent linear distribution will be continuous. Determining the lift 
distributions for the cascades is essentially a stepwise process in that 
the vortex drag criterion is satisfied over part of the space and then 
the wave drag criterion over part, alternating back and forth until both 
conditions are satisfied everywhere. In this example rotational symmetry 
is assumed and the center cascade is used to satisfy the vortex drag 


re quirement s ; thus , 


the outer region —= < r < R of this cascade must 

ITS 
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carry a constant intensity of lift. The cascades of radius R/lf2" are 
used to give the equivalent linear distribution the required elliptic 
shape for R/i/ 2 - ^ £ RV2". The next step is to evaluate the distribution 

over another section of the center cascade to give constant lift inten- 
sity when elements are summed up in the free stream direction, then 
satisfy the wave drag criterion with the next cascade, etc. This proc- 
ess is continued working inward to the center of the space; although 
an infinite number of cascades are required the total wing area is 
finite. Each of the small cascades has a’ radius l/)f2 times the radius 

of the next larger one and the total wing_area is 3 = 2.172itR 2 (Ch. VI B) 
It should be noted that this is not necessarily the minimum wing area 
that could be used, so the distribution obtained is an optimum one with 
respect to wave and vortex drag only and not with respect to friction 
drag. 


B_j_ THE OPTIMUM DISTRIBUTION OF LIFT THROUGH AN ELLIPSOIDAL SPACE 

The spherical space with Its optimum lift distribution can be 
changed into an ellipsoidal space with a corresponding lift distribution 
by a scale transformation of one of the cartesian coordinates. This 
transformation transforms planes into planes so that elliptical loadings 
are preserved for the ellipsoid and minimum wave drag; is obtained. 

Also a constant intensity of lift over the wake cross-section is 
maintained for the ellipsoid so that the vortex drag_is also a minimum. 

Although the optimum lift distribution for an ellipsoid is obtain- 
able from the spherical case, the value of the minimum drag Is not nec- 
essarily the same. For an ellipsoid formed by revolving an ellipse of 
semi -major axis B and semi -minor axis R about the free stream (major) 
axis, the optimum distribution of lift is • 


''opt 


1 - (x/B) 2 - (y/R) 2 - (z/r) 2 ! 


1/2 


The wave drag, computed by Hayes’ method, is 
*Wn wave 


P 2 L 2 


8*qR 2 [(B/R) 2 + p 2 J 
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and the vortex drag is also a minimum. 


Dmin vortex = 



8k qR^ 


so that the total drag is 


•Pmin - 


8k qR^ 


J -P 2 

[(b/r) 2 + p 2 ] 


+ U 


For B = R the results reduce to the spherical case . 


Several limiting cases can he examined; in one an ellipsoid is 
collapsed into a horizontal planar wing of elliptic planform carrying 
constant pressure. Optimum cases of this type were first discussed hy 

R. T. Jones Another limiting case which gives minimum drag occurs 
when an ellipsoid is collapsed into a plane normal to the flow direc- 
tion (b/r — > 0). Then the wing system can he interpreted as a uniformly 
loaded airfoil cascade (of zero chord and gap) within the elliptical 
cross-section. The entire cascade can he analyzed as a two-dimensional 
system. If the chord is chosen to be p times the gap then the airfoils 
in the cascade are non- interfering hut the lift distribution is suffi- 
ciently continuous (Fig. 6b-l) . In other words, when the cascade is 
cut hy planes Inclined at the Mach angle, the resulting load distribu- 
tions used in Hayes 1 method will he continuous. The total wing area is 
then (3 times the area of the ellipse.. 
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Fig. 6b-l: Examples of airfoil spacing in cascades 
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A third limiting case is the slender "body obtained when B/R — > »; 

then 


2 22 - 

^min = ^ p + o = ^Vortex + ^wave 

2*q(2R) 2*q(2B) 1 

The wave drag portion is the same as that obtained by Jones for a planar 

slender wing^ 1 -^, while the vortex drag for. the spatial distribution is 
one -half that obtained by Jones for the planar distribution. 


C . THE OPTIMUM DISTRIBUTION OF LIFT THROUGH A "DOUBLE MACH CONE" 


Consider a space consisting of two Mach cones placed base to base 
(Fig. 6c-l) . If a uniformly loaded cascade of airfoils (with zero gap 
and chord) is placed at the maximum cross -section of this space then 



Fig. 6c-l: Double Mach cone space 

with optimum cascade 


elliptic loadings will be obtained when the space Is cut by planes 
inclined at the Mach angle. This airfoil cascade consequently produces 
the minimum possible wave drag for wing systems contained within the 
space and carrying a specified lift. The uniform distribution of load 
over the circular cross-section Insures minimum vortex^ drag also, so 
the lift distribution is an optimum for the double Mach cone. 

The value of the minimum wave drag (obtained by Hayes' method) Is 
D wave = i — — — ^ and the vortex drag has the same magnitude in this 

d Jtq(2Rr 


case . 
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The wave plus vortex drag is then D = — . This is equal to 

«q(2R) 2 

the minimum vortex drag alone for a -planar wing of span 2R. If the air- 
foil cascade is compared to the largest planar wing of diamond planform 
which can he contained within the. double Mach cone, the minimum wave 
plus vortex drag of the diamond planform is approximately 1.52 times 

greater than for the cascade^). 

Again It must he emphasized that the drag comparison is not com- 
plete without the inclusion of viscous drag and thickness drag for the 
wing system. 

Since the circular cascade is an optimum arrangement, it satisfies 
Jones' criterion (Ch. V). This can he checked as follows: By two- 

dimensional analysis the downwash, €, in the aft Mach cone is 2a where a 
is the angle of attack of each airfoil (Pig. 6c-2) . Since the downwash 
is zero in the fore Mach cone, the downwash velocity in the combined 
field Is constant and equal to aU throughout the double Mach cone. 



Fig. 6c-2: Two-dimensional analysis of downwash in rear Mach cone 

of an optimum cascade 
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Far behind the cascade in the wake of the wing system s = a; this 
can be shown by equating lift to rate of change of vertical momentum. 

The individual wings of the cascade are non-interfering and, in the 
limit as gap and chord go to zero, have two-dimensional wing character- 
istics. The wing area for a sufficiently continuous lift distribution 
(Ch. VIB) is equal to the cascade cross-sectional area A times (3.' Con- 
sequently L = C^qS = (4<x/|3)q(|3A) . By Munk's criterion (see Ch. V and 

Ref. 12) the downwash in the Trefftz plane over the area behind the 
cascade is constant; thus, the vertical momentum, of the fluid in the 
downwash region behind the cascade is (pAU) (eU). The vertical, momen- 
tum of the surrounding fluid can be evaluated from the known "virtual 
mass" of a solid circular cylinder of cross-sectional area A moving 
downward in the fluid; this latter momentum is equal to that of the 
downwash region itself. Thus, by the momentum theorem, L = 2pAU(eU) 
and equating the two expressions for L gives e = a. 

The airfoil cascade is not the only distribution of lift in the 
double Mach cone which has minimum wave drag. A true lineal distribu- 
tion of lift distributed as an elliptic loading along the axis of the 
double Mach cone will produce the same minimum value of wave drag. So 
also will a lift distribution of constant intensity throughout the entire 
double Mach cone. However, the latter two cases will not give the mini- 
mum value of vortex drag; in fact, the true lineal distribution will 
have infinite vortex drag. 



A sufficient condition for minimum drag is that each equivalent 
lineal distribution of lift should be elliptic (Ch. V). For the spherical 
space the Be equivalent lineal distributions will be the same at a.1 1 angu- 
lar stations if the optimum lift distribution is rotationally symmetric. 
For simplicity, examine the problem from the angular position 9 (on the 
control surface) equal to 90°; then the Mach planes will be parallel to 
the yaxis. The notation to be used is illustrated in Fig. a6-1j cylin- 
drical coordinates (£, S,0) and the radial coordinate r will be used. 

If the spatial lift distribution is l( r) = z(^£ 2 + S 2 ) then the equiva- 
lent lineal distribution along the £ axis will be 


>2 ,2 


-r r <^) 

J s=o J rf=n ' ' 


t 2 + S 2 )s d0 dS = 2n 




Sl[H 2 + S 2 )dS 


MACH PLANES FOP G = 



?+rf=R a -S* 


Fig. A6-1 


However, 


= Kl/l - (t/R ) 2 
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where K depends on the total lift of the sphere . Introducing the radial 
coordinate r, the integral equation to be solved is 


Kjjl - (£/R) 2 = 2jt J rl(r)dr 


The solution to this equation, found by differentiation with respect 
to £, is 


l(r) 


K 

2jtR^R 2 - r 2 


The total lift of the sphere is 

L = 2 r k|/i - (C/R) 2 d£ = 

so that the distribution of lift for mi nimum wave drag is 




For application of Hayes ' method, the equivalent lineal distribution 
along the x axis is needed. A plane £ = £ ' intersects the x axis at 
x = -M£ ' ; since the distribution is spread- out over a larger distance 
along the x axis, its maximum intensity will be less; thus. 


F < x )opt = If 1 - (x/MR) 2 = |jy|l - (x/MR) 2 


Hayes defines two functions such that for the lifting case (Ch. IV) 


F = 


P U 6 Z 

P 


h = -g z £in 0 


-2|3L sin 0 
itpUMR 


1 - (x/MR) 2 
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The expression for the wave drag contribution at each angular station 0 
is, from Eq.. (4-C-30), 


dD 

d0 


-P 

8* 2 


JJ h '( x 2) h '( X l) Zri | X 2 - x-^jdx^ dx £ 


and the total wave drag is 



de 

d0 


The integration for dD/cL9 has been carried out by Sears in terms of 
a Fourier series expansion of an arbitrary function h. For the wave 
drag optimum the distribution h is elliptic and only the first term in 
the series for h appears. (Note the similarity to the vortex drag opti- 
mums in incompressible flow.) If h = cfl - (x/MR) 2 then dD/d0 = pC 2 /l6. 
Substituting in the equations above leads to the final result, 

r _ 0 21 * 2 
^wave - ~z 

8jtqR4^ 
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CHAPTER VII. THE OPTIMUM DISTRIBUTION OF VOLUME ELEMENTS ALONE* 


Aj_ THE SINGULARITY REPRESENTING AN ELEMENT OF VOLUME 

The investigation of lift distributions is simplified by the use of 
a singularity which represents an element of lift. This singularity is 
the elementary horseshoe vortex. The intensity of lift corresponds to 
the strength of the singularity and the location of the lift force is 
identical with that of the bound vortex. The study of volume (or 
thickness) distributions is similarly simplified by identifying the sin- 
gularity which corresponds to an element of volume . 

Consider a source and sink of equal strength and located on the same 
streamwise' line. In each unit of time a certain quantity of fluid is 
introduced into the flow pattern by the source and the same quantity is 
removed by the sink. The volume occupied by the fluid flowing from source 
to sink depends on the strength of the source and sink and the distance 
between them, and also depends on the velocity and (tensity of the fluid 
flawing from source to sink. However, if the volume is to be considered 
a linear function of the strength of the singularities, then the mean 
value of density times velocity must be unaffected by the perturbation 
velocities created by the source and sink. This means that in a line- 
arized treatment of the problem the fluid flowing from source to sink 
may be considered to have free stream density and velocity. 

Let m = Mass of fluid introduced per unit time 
d = Distance between source and sink 
p Q = Free stream density 

U Q = Free stream velocity 

Then the volume occupied by the fluid is 
vol = md^p 0 U 0 ) 

Since the volume is proportional to md, doubling the intensity of 
source and sink and halving the distance between them should produce a 
shorter, but thicker volume of the same magnitude. This suggests pro- 
ceeding to the limiting case (as in incompressible flow) where the source 


* 


The contents of this chapter have appeared in the paper "The Drag of 
Non-Planar Thickness Distributions in Supersonic Flow," published in 
the Aeronautical Quarterly, Vol. VI, May 1955* 
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and sink are combined in a dipole with, axis in the* free stream direction. 
This singularity should represent an element of volume, although the 
fineness ratio of the element is zero. 

The potential for a unit source at (|,0) in supersonic flow is 

rf>Q = ================ 

2s||(xT|)2T 


where fS = |fM 2 - 1; x and | are coordinates in the streamwise direction 
and r is radial distance from the x axis. 

Differentiating with respect to x gives 


X_ 17 v2 o 2-13/2 

2jt | (x - |) 2 - prr 2 j 

where is the potential for the unit dipole or an element of volume 
equal to l/U Q . 

B. THE DISTRIBUTION OF VOLUME ELEMENTS 


For a distribution of volume elements along the | axis with inten- 
sity f(|), starting at | = 0, the potential is 

a m i_ r x ' pr f(s)U - i)de 

251 Jo r . 2 2 21 

|_(x -ft ) - ^r 2 J 

Integration by parts gives 


f(!) 

x- S r L r x- ? r f , (5)1| 

f(x - 1)2 - P 2 r 2l 

0 2lt Jo ( l(x - 1)2 - p2 r 2 



The first term in the expression for the potential is infinite, 
and apparently corresponds to the "roughness" of the body, which iB an 
assembly of blunt elements (see illustration) . 


a 
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The smoothly faired body (indicated by dash lines) is all that we 
are concerned with, and this creates the finite part of the potential. 

This finite part is also the potential for a source distribution of 

intensity equal to +f'(j-). This source distribution can be used to con- 
struct a body of revolution extending from -l / 2 to +l/2. 


The shape of the body of revolution created by the singularity dis- 
tribution may be obtained approximately by slender body theory or more 
accurately by "exact" linear theory. In the first case the volume is 

f (e 

. -15J—5 - which agrees exactly with the sum of the volume elements. 


/ 


- 1/2 


U r 


An example of the second case is shown in Fig. 4f -2 where a singularity 
distribution on the axis is interpreted first by slender body theory 
then by "exact" linear theory as a "bump" on a cylinder. The bump shapes 
and locations are quite different but the volumes are identical. This 
has been proved by Lagerstrom and Bleviss and generalized by Bleviss in 
Ref. 22. ' 


A planar distribution of volume elements may be interpreted by 
("exact") linear theory as a thin planar wing. The volume contained in 
this wing is exactly equal to the sum of the volume elements. 

The concept of the volume element is not necessary for the study of 
smooth slender bodies of revolution and planar wings, since these con- 
figurations are relatively simple . However the use of the volume element 
does help to clarify problems involving more general spatial distributions 
of thickness. — - 

The points to be emphasized are that fixing the sum of the volume_ 
elements fixes the total volume, and fixing the distribution of volume 
elements determines the drag. It is therefore possible to study the 
drag of a distribution of volume elements without calculating the exact 
shape of the corresponding body. This is analogous to the fact that the 
drag of a distribution of lifting elements can be studied without calcu- 
lating the twist and camber of the corresponding wing surfaces. 


C. THE DRAG OF VOLUME DISTRIBUTIONS ON A 

STREAMWISE LINE AND THE SEARS-HAACK BODY 


A body of revolution may be constructed from a distribution of vol- 
ume elements along a streamwise line, or from the equivalent distribu- 
tion of sources. The body constructed from volume elements Is an 
"Infinitely rough" body and has infinite drag. However, discarding the 
infinite part of the potential leaves a "smooth" body (with finite drag) 
which is equivalent In every respect to the body created by a source 
distribution. 
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If f (x) is the intensity of tie volume element distribution for a 
body of revolution of length "Z" then the drag is given by^^ 

r> + Z /2 p + Z /2 . 

D - £j_ l/2 J_ i/2 f "( X l) f "( X 2) !ll K - ^2 


To maintain constant total volume according to linearized theory 
p +Z/2 

it is necessary that / f(x)dx = Constant. The body shape giving 

J -Z/2 

minimu m drag for a given length and volume has been determined by Sears 

and Haack^ 1 ?) independently. The corresponding f(x) (which is propor- 
tional to the cross-sectional area) is given by 


(16) 


16 

l . /^ 2 1 

1 3/2 J 

p + J /2 

1 f (x)dx , 

-Z/2 _ 8LTo volume 

'i - 

3* 

\i) J 

z 3 * z /2 

UJ J 


3/2 


Thus the optimum distribution of volume elements along the axis 
corresponds to the cube of an elliptical distribution. (For lifting 
elements the optimum distribution is elliptical.) 

The value of the minimum drag is 


%in = & 


o+Z/2 

2 


/ f(x)dx 


r “i 

J -Z/2 

8q/lf 

volume 

U 0 (z/2) 5 

« w 

_(Z/2) 5 _ 


D. THE SEARS-HAACK BODY AS AN OPTIMUM 
VOLUME DISTRIBUTION IN SPACE 


If the volume elements are not confined to a single streamwise line, 
then the drag contributions at different angles, 0, on Hayes 1 .cylindrical 
control surface are not necessarily the same. For any one angle, 0, the 
drag is given by 
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Here f(x,0) is determined by the use of "Mach planes" for the 
angle 9. All the volume elements intercepted "by any one "Mach plane" 
are transferred (in the plane) to the streamwise axis. The resulting 
distribution along the axis is f (x, 9 ) . The problem of finding the mini- 
mum drag contribution at the one angle £ is then similar to the Sears- 
Haack problem. If f (x, 9 ) corresponds to the cube of an elliptical dis- 
tribution for every 9, then the total drag is a min imum , and the drag 
contribution at each 9 is a minimum and corresponds to that of an equiv- 
alent Sears-Haack body. 

It Is not always possible to simultaneously minimize the drag 
contributions at all angles 9 . However If we consider the optimum 
distribution of thickness 
within a space which has 
rotational symmetry about 
a streamwise axis, then 
it may be possible that 
all the equivalent bodies 
are Sears-Haack bodies 
having the same length. 

For example, consider that 
a double Mach cone bounds y 

the space within which < 

thickness is to be distrib- 
uted. The Sears-Haack 
body placed on the axis 
is an optimum for this 
space . It has the same 

drag contribution at SEARS-HAACK body bounded 

every angle on the cylin- BY double suach cone Space 

drical control surface, 
and of course, the 

"equivalent" body of _ - 

revolution for any angle 9 is identical with the real body. However, 
a "ring" wing (which carried no radial forces) plus a central body of 
revolution can be designed 
to have exactly the same 
drag as the Sears-Haack 
body. The equivalent 
bodies of revolution are 
all identical with the 
Sears-Haack body. This 
is discussed in the next 
section. (For the case 
in which radial forces 
are carried on the ring 
wing see Ch. IX.) 



BODY HA V/N<3 SAME DR A <3 
AS SEARS-HAACK BODY 
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E. RING WING AND CENTRAL BODY OF REVOLUTION 

COMBINATION HAVING THE SAME DRAG AS A SEARS-HAACK BODY 

Consider a ring-wing plus a central body of revolution contained 
within the space bounded by a double Mach cone. Because of the rota- 
tional symmetry of this particular system, the equivalent body of revo- 
lution is independent of the angle 0 on the cylindrical control surface. 
In this case, if the local radial force on the wing is everywhere zero, 
the drag of the equivalent body of revolution is, according to Hayes' 
formula, identical to the drag of the original system. Thus, a ring- 
wing (which carries no radial force) plus a central body of revolution 
will have exactly the same drag as a Sears -Haack body if the equivalent 
body of revolution is a Sears -Haack body. 

To design such a system, we may select any smooth, slender profile 
for the ring-wing and compute the cross-sectional areas cut from this 
wing by a set of parallel Mach planes. These areas must then be sub- 
tracted from the cross-sectional areas which would be cut from a central 
Sears -Haack body by the corresponding Mach planes. The resulting area 
difference defines the area distribution (in tbs Mach planes) of the 
correct central, body. (This area must be projected normal to the flow 
direction to obtain the cross-sectional area of the central body defined 
in the us ual way. ) This body, together with the ring -wing originally 
selected, is an optimum distribution of thickness within the double Mach 
cone space. 

As an example, consider a ring-wing with thickness distribution 
corresponding to a bi-parabolic arc profile. The camber necessary for 
zero local radial force need not be determined, since it does not affect 
the shape of the central body. Assume that the wing is six percent thick 
located half-way between the axis and the apex of the space. If the 
central body of revolution is designed so that the equivalent Sears-Haack 
body is of fineness ratio 5, the resulting shape of the central body of 
revolution is as shown in Fig. 7e-l. 
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Fig. 7e-l: Cross-sectional view of ring-wing and central "body 

(an optimal distribution of thickness within the double Mach 
cone space) 

F\ OPTIMUM THICKNESS DISTRIBUTION FOR A 
PLANAR WING OF ELLIPTICAL PLANFOEM 


It is desired to find the optimum thickness distribution for a 
planar wing of elliptic planform and given volume; this problem was 

first solved by R. T. Jones^ 1 ^). A geometrically simpler problem, which 
will be examined first, is to find the optimum thickness distribution 

for a circular wing of given volume. The method of Hayes ^ in which 
the drag is evaluated by summing increments of drag at each angular sta- 
tion around a cylindrical control surface far away from the body, will 
be used. For the total drag to be a minimum, the increment of drag at 
each angular station should also be a minimum. 

If the thickness distribution of the circular planform is rota- 
t tonally symmetric, then the equivalent bodies at each angular station 
will have the same shape (although different "fineness ratios") due to 
symmetry. If t(r) is the thickness distribution to be optimized for a 
given volume V, then __ 
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V = 2rt 



t(r)r dr 


(7f-l) 


where R Is the radius of the circular wing and r, 0 are polar coordinates 
from the wing center (Fig. 7f-l). The area cut out at each point along 

the 1 axis by planes normal 
to that axis is 



The equivalent lineal distribution along the x axis is 


r R 

S(x) = 2 cos p / 

J v r 


rt(r )dr 


x cos p jfr 2 - x 2 cos 2 p 


(7f-5) 


with 


'+R sec p 


/ ' S(x)dx = V 

“ -R sec p 

For minimum drag, this distribution should be (Ch. VF) 

r 2 i 5/a 

S(x) a 11 - (x/R sec p) 

Thus the integral equation to be solved for t(r) is 


k|i - (x/R sec p) 2 ^J = 2 cos p J 


t(r)r dr 


x cos p x 2 cos^p 


(7f-4) 


(7f-5) 
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where K is a constant dependent upon the given wing volume . By a suit- 
able transformation of coordinates, Eq. (7^-5) ms-y written in the form. 


K ^ t (a)dq 

H sec V- ~ J 0 ||0 - a 


(7f-6) 


where 


0 = 1 


(x/R sec n) 


2 


a = 1 - (r/R) 2 


Eq. (7f-6) is called Abel's equation and its solution is well known, 
c.f.. Ref. 19. The solution to Eq. (7f-6) is 


t(r) 


3K 

iffi cos \l 


1 - (r/R) 2 _ 


and substitution of this in Eq. (7f-l) determines K; then 


. t(r ) 



(TC-7) 


Equation (jf-7) thus gives the distribution of thickness which will 
result in minimum drag for the circular planform wing of given volume,. 

To apply the circular planform solution to the original problem of 
finding the optimum thickness for an unyawed elliptic planform, make 
the following change of coordinates: _ 



Fig. 7f-2 
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The circular wing is then transformed into an elliptic wing whose equa- 
tion is 



It can be verified that the thickness distribution 


t 



(if -9) 


obtained frcsn Eq. (7f-7) through the transformation Eq. (7f-8) is the 
optimum for this more general case; that is, the equivalent linear dis- 
tribution for Eq. (7f-9) with a set of Mach planes inclined at the 
angle n as shown in Fig. Jf-2. is 


where 



(7f-10) 


l = |/a^ + b^tan^p. 

Since Eq. (7f-10) represents a Sears-Haack body, the thickness given by 
Eq. (7f-9) is optimum for the unyawed elliptic wing. 

Determination of the total drag in this optimum case involves an 
integration of the drag increments from these Sears-Haack bodies as 
seen at each angular reference station. If the reference station is 
at an angle 0 from the horizontal, then the Mach planes cut the elliptic 
planform at an angle p defined as (Ch. IVC). 


tan |j. 



1 cos 0 


(7f-H) 


D = 



and the total drag is 
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The increment of drag at each reference station is (Ch. IVC) 


dD _ _ _q_ 

d9 kiiu q u o 

4qV 2 


J S ,, (x)S"(S)ln|x - || dx dS 


(TC-12) 


* 2 Z^ 


and the total drag for the optimum thickness distribution Eq. (7f-9) is 


Defining 



and 


D - C-Qqitab 


then 



This result agrees with that 


given by Jones 


(14) . 


(Tf-15) 


(7f-l4) 
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CHAPTER VIII. UNIQUENESS PROBLEMS FOR OPTIMUM DISTRIBUTIONS IN SPACE 


A^ THE NON -UNIQUENESS OF OPTIMUM DISTRIBUTIONS IN 
SPACE - ’'ZERO LOADINGS 11 


In the subsonic flow of a perfect fluid the only drag caused by a 
lifting wing is vortex drag. The minimum possible vortex drag for a 
planar wing is obtained when the spanwise lift distribution is elliptical. 

According to Munk's stagger theorem/^) the chordwise location of the 
lifting elements is unimportant, so there are infinitely many distribu- 
tions of lift over a given planform which produce the minimum drag. 

In supersonic flow lift causes both vortex drag and wave drag. The 
chordwise location of lifting elements is still unimportant in deter- 
mining vortex drag, but does affect the wave drag. For this reason the 
optimum lift distribution for a planar wing is generally unique in super- 
sonic flow. However, spatial lift distributions offer more freedom in 
the arrangement of lifting elements and the optimum distributions in 
space are not generally unique even in supersonic flow. 

For example, the minimum wave drag due to lift in a double Mach 
cone space can be attained with each of three different simple lift dis- 
tributions. (See VI-C.) The first is a constant intensity over the 
circular disc located at the maximum cross-section of the space. The 
second is an elliptical intensity concentrated on the axis of the double 
Mach cone. The third is a constant intensity throughout the entire 
double Mach cone . If the first two distributions are superimposed, one 
carrying a unit of positive lift and the other a unit of negative lift, 
the result is a net lift equal to zero. Also, the net strength of the 
lifting elements intercepted by any cutting plane inclined at the Mach 
angle is zero. This means that the combined distribution has zero wave 
drag. Furthermore, there are no disturbances whatsoever produced on 
the distant control surface near the Mach cone and no wave drag inter- 
ference can exist with any other loading. If another such combined dis- 
tribution with opposite sign is placed on the same streamwise line with 
the first one, then, by Munk's stagger theorem, the vortex drag is zero 
also. This is one example of a "zero loading" (see illu stration), and 
many others can be constructed. 
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A "ZERO LOADING” PLACED W/TR/N A/V ELL/PSO/DAL SPACE 

Such a "zero loading" placed within any space alters neither the 
lift nor the drag of the orig i nal lift distribution. For this reason 
optimum lift distributions in three dimensions are never unique (unless 
the space degenerates into a surface ) . 

Similar arguments can lie applied to optimum distributions of volume. 
For an example of. non-uni queness in such cases see Ch. Vll. 


B. UNIQUENESS OF THE DISTANT FLOW FIELD ' " _ 

PRODUCED BY AN OPTIMUM FAMILY 

It has been shown, that optimum lift or volume distributions in 
space are not generally unique, since a group of optimum distributions 
can be obtained from one given optimum distribution by superposition of 
"zero loadings." Each member of the group produces the - same (minimum) 
value of drag for a given total, lift or volume . ~ 

From the method of construction of this group (by the use of "zero 
loadings") it follows that each member produces the same velocity per- 
turbation field in the Trefftz plane and on the distant control surface 
near the Mach cone. It can also be shown that there are no optimum dis- 
tributions outside this group, since all possible optimum distributions' 
are indistinguishable from the "distant" viewpoint. 

Assume that f l opt £) and - f^pt^' 1 ^^ are members of the opti- 
mum family not included in the original group (whose members were related 
through "zero loadings"). Assume also that and fp 0 pt do not P ro_ 

duce identical perturbation velocity fields far from the singularity “ 
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diet ribut ion. 


Then the drag of f-< equals the drag of f 2 ( 

x opt ^opt \ 


or 


^lopt = ^2opt ) ^ definition of the optimum family. Also fp^^ ma ^' 

set equal to f n . + AC. where Af carries zero net lift (or volume), 
J -opt 7 

hut has a velocity perturbation field which is not identically zero far 

from the singularities. 


The distribution Af is orthogonal to (does not interfere with) fp Q ^. 

This follows because any given lift or volume distribution can be improved 
through combining it with a distribution having zero net lift or volume 
if there is interference drag. However f by definition, cannot be 

improved, and must, therefore, be orthogonal to Af . 


D 2opt = D lopt + D Af> but we also 

know that Do = D-i and. therefore. D A -p must equal zero. Here we 
*~opt x opt * ^ 

can obtain a contradiction since both the vortex drag and the wave drag 
depend on the squares , of velocity perturbations (in the Trefftz plane 
and far out on the Mach cone) and the drag contribution from each portion 
of the control surface is non-negative. If Af produces any disturbances 
far from the lifting system it must have positive drag, and so Af must 
produce Identically zero disturbances to have zero drag. 

The above contradiction shows that all the members of the optimum 
family are indistinguishable from the distant viewpoint. 


Since Af L is orthogonal to f 


lopt* 


If drag Is computed from the "close" viewpoint the above argument 
cannot he made. Drag contributions then appear as the product of local 
pressure times angle of attack on the wing surfaces, and these quantities 
are not necessarily non-negative at every point on the surface. 


UNIQUENESS OF THE ENTIRE "EXTERNAL" FLOW FIELD 
PRODUCED BY AN OPTIMUM FAMILY 


It has "been shown that any two members of an optimum family produce 
Identical velocity perturbations on the distant control surface. 


If f-L (6,11,0 and f 2 (6,11,0 are two members of an optimum 
opt opt 


family, then f 


\>pt 


-opt 


must produce identically zero velocity per- 


turbations on the distant control surface, and the drag will be zero. 


Let "S" designate the space within which the singularity distribu- 
tion f -i - fp exists, and let "E" represent the external flow field 
- L opt c opt 
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consisting of points whose aft Mach cones 
do not intersect "S." Assume that at 
some point in the external field "E" the 
resultant velocity vector is inclined to 
the free stream direction. Then an ele- 
mentary wing can he inserted at that 
point with the angle of attack adjusted 
to give negative drag on the wing. Since 
the singularities in "S" are outside the 
aft Mach cones of all points in "E, " the 
net drag change produced by the elementary 

wing is negative. However, f-, - f 0 

x opt ^opt 

is a singularity distribution causing zero 
drag, so f lopt - f 2opt P lus the elementary 

wing is a system having negative drag, although it is an isolated system 
inserted in a uniform flow field. However, the drag of this system eval- 
uated on a distant control surface ccmes from a summation of positive 
quantities and cannot be negative. This contradiction ahows that the 

external flow field "E" produced by f-, - fo must 'consist of velocity 

x opt c opt ~ 

vectors aligned with the free stream direction. These vectors must also 
have the magnitude of the free stream velocity; hence, the external flow 
field is completely undisturbed, and it can be concluded that al 1 members 
of the optimum family produce the same flow pattern in the external 
field "E." 

It is of interest that a similar proof cannot be made for subsonic., 
flows . In such cases there is no external region where an elementary 
airfoil can be inserted without producing interference effects at the 
original singularities. • - 



II EXISTENCE OF SYMMETRICAL OPTIMUM DISTHIBbTIONS ~ 

IN SYMMETRICAL SPACES ' .H. 1" 

It can be shown that, if the boundary of a space has a horizontal 
plane of symmetry, then there is one member of the family of optimum 
lift distribution within the space which is symmetrical, about the plane. 
The proof is as follows : “ 

Let * 0 pt(x,y,z) represent an optimum lift distribution in the space. 
The distribution Z 0 pt( x >y> -z) has the same drag and lift" (the drag of 



NACA TM 1421 


the individual lifting elements 
is unaltered hy the change of 
position, and the interference 
drag of any element pair is 
unaltered also). 

Since (x,y, -z) has the 
same lift and drag as l Q ^(x } y,z) 

it is also a member of the optimum 
family. All members of the opti- 
mum family produce the same exter- 
nal flow field, and any distribu- 
tion producing that field is an 
optimum. The distribution, 

7^opt(x,y,-z) + |l 0 pt(x ; y ; z) pro- 

duces the same external flow field 
as Z 0 p*t(x,y, z) . It is, therefore, an optimum, and since it is also 
symmetrical about the horizontal plane the proof is completed. 

Similar proofs can be developed for cases where lift, thickness 
and side force elements are present. Also certain other planes of 
symmetry can be used. 
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CHAPTER IX. INVESTIGATION OF SEPARABILITY OF LIFT, 
THICKNESS AND SIDEFCRCE PROBLEMS* 


A. THE SEPARABILITY OF . OPTIMUM 

DISTRIBUTIONS PROVIDING LIET AND VOLUME 

Separability Questions 

For the purpose of drag evaluation a complete aircraft is repre- 
sented by a distribution of lift - elements, volume elements and possibly 
sideforce elements in space. A certain net lift must be provided to 
support the weight and a net volume must be provided to house payload, 
fuel, structure, etc. The drag should then be made' as small as possible 
with the net lift and volume equal to the prescribed, values. 

Several questions arise . Can we first study the problem of how 
best to provide the required lift (with no net volume), then determine 
the best way to provide the required volume (with no net lift), and 
finally by superposition obtain the optimum distributions of singulari- 
ties for simultaneously providing the net lift and volume? If this 
procedure is possible will the drag of the combination be the sum of 
the drags of the two superimposed distributions ? Does the optimum way 
of providing the lift with no net volume require only lifting elements 
or are volume and sideforce elements necessary? Similarly does the 
optimum way of providing the volume with no net lift require singulari- 
ties other than volume elements? 

For horizontal planar systems the answers to these questions are 
comparatively simple. The lift and volume problems can be studied sepa- 
rately and the optimum singularity distributions superimposed. The drag 
of the combination JLs the sum of the drags of the individual distribu- 
tions. Finally, the optimum way of providing the lift requires only 
lifting elements and the optimum way of providing volume requires only 
volume elements . 

All of the above results follow from the fact that in horizontal 
planar systems there is no interference drag among lift, sideforce, and 
volume elements. However this is not true in general for non-planar 
systems, and consequently the above problems must be re-investigated 
for these more general configurations. 


* Portions of this chapter have appeared in the paper "The Drag of 
Non-Planar Thickness Distributions in Supersonic Flow," published in 
the Aeronautical Quarterly, Vol. VI, May 1955- 
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Optimum Distributions Providing Lift and Volume 

In non-planar distributions of lift, sideforce and thickness there 
is generally interference among the different singularities. This means 
that the drag for a given net lift may in some cases he decreased by 
adding thickness or sideforce elements and taking advantage of negative 
interference drag. 

In order to study such cases let 7,(x,y,z), t(x,y, z), and s(x,y,z) 
represent respectively distributions of lift, thickness and sideforce 
in x,y, z space within some boundary. Here we will exclude, without loss 
of generality, those distributions of l and s which are completely equiva- 
lent to elements of volume or thickness (see III -A, the closed vortex 

line). Let z) + t-j_ 0 (x,y,z) + s-j_°(x,y,z) give the minimum pos- 

sible drag for one unit of net lift and zero net thickness and sideforce. 
(The superscript simply indicates the net lift or thickness or sideforce 

of the distribution.) Also let l 2 °(x,y,z) + t 2 ^(x,y,z) + s 2 °(x,y,z) give 

the minimum possible drag for one unit of net volume and zero net lift 
and side force. We ask what distribution gives the minimum drag when 
both the net lift and net volume are simultaneously prescribed and equal 
to L 0 and V Q respectively? 


Consider the distribution A(x,y, z 



+ tn 



which gives the prescribed net lift and volume. For this to be the opti- 
mum it is necessary and sufficient that it be orthogonal to every distri- 
bution 1° + t° + s°, which contains zero net lift, zero net volume and 
zero net sideforce. For example, any loading 1° + t° + s° multiplied by 
an arbitrary constant C can be superimposed on A without altering the 
net lift, L Q and net volume, V Q . If this distribution 1° + t° 4* s° were 

not orthogonal to A, then C could be adjusted to give a negative inter- 
ference drag with A greater than the drag of C(Z° + t° + s°) by itself. 
Hence the distribution A could be improved and therefore would not be 
an optimum. It is also true that any possible improvement of A must be 
obtainable by superposition of a loading of the type C(l° + t° + s°) on A. 
So for A to be an optimum it is both necessary and sufficient that A be 
orthogonal to any loading 1° + t° + s°. 


However, + ^i° + s i° ^2° + *^2^ + s 2° are eac k orthogonal to 
any 1° + t° + s° since each one is an optimal distribution in its own 
restricted class. Therefore because of the linearity of the interference 

terms + t^° + s l°) + ^o(^2° + ^ 2 "*" + s 2 0 ) is orthogonal to any 

1° + t° + s° and A(x,y,z) is the optimum distribution having lift = Lq 
and volume = V G . 
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The Drag of tbs Optimum Distribution Providing Lift and 
Volume in a Region Having a Horizontal Plane of Symmetry 

The drag of the optimum distribution A(x,y, z) must next be deter- 
mined. If Zj 1 + t-j_° + S 2 _° is orthogonal, to Zg° + tig 1 + s 2 ° then the 

drag of A(x,y, z) is just the sum of the drags of Lq (z^ 1 + t^ 0 + s^ 0 ) 
and of V 0 (z 2 ° + + s 2 °). We knov that (z 2 ° + tg 1 + s 2 °) X (b^ 0 + s-^ 

so the question arises is (zg° + + S2°) also J-Zi 1 ? (Here the sym- 

bol 'X' indicates orthogonality.) 


In order to answer this question it is convenient to represent 
l^(x,y,z) by a concentrated lift of one unit Z^g 1 plus a distribution 
containing zero net lift Z-^°(x,y, z). The concentrated unit of lift can 
be placed anywhere in the space and then Z 2 _°(x,y, z) is simply the dif- 
ference between Z^^y, z) and Z^g 1 . Similarly it is convenient to 
replace t^g by t^gg + t° 2 (x,y, z). The optimum distribution is then 


A(x,y, z) = L 0 (z^g 1 + Z X ° + t-L 0 + s-^ + V Q (tgg 1 + Zg° + tg° + Sg°) 
The distributions in brackets are orthogonal if. 


Z*j^g _1_ tgg + Zg 4* tg + .Sg^* 


or if 


tgg 1 X Zqg 1 + Zj_° + t ± 0 + Bq° 

The concentrated unit of lift Z 1 g 1 may be located at any point in 

the space and has the same interference drag with ("tgg 1 + Zg° + tg° + s 2 ° 
for all locations . Thus if there is any point in the space where a unit 
of lift has no interference with (tgg 1 + Zg° + tg° + s 2 0 ^ the orthogo- 
nality of the two components of A(x,y,z) is assured. (This does not 
depend on the connectivity or the convexity of the space.) 
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For example, if the boundary of a space has a horizontal plane of 
symmetry, then there are optimum distributions in the space having sym- 
metry properties. (The proof is similar to that given in Ch. VIII for 
a lift distribution.) If some portion of the plane of symmetry is con- 
tained inside the space then the concentrated unit of lift can be located 
in this plane and orthogonality demonstrated. 


B. THE NON - BlTKHi’KKENCE OF SOURCES WITH OPTIMUM 

DISTRIBUTIONS OF LIFTING ELEMENTS IN A SPHERICAL SPACE 


In general there is interference between non-planar distributions 

of sources and lifting elements, as shown by Hayes . This means that 
in general the optimum distribution of singularities which provides one 
unit of lift may contain volume elements or sources as well as lifting 
elements . However for certain spaces it can be proved that there is no 
interference between a source and the optimum distribution of lifting 
elements alone. So for these spaces the optimum way of providing lift 
requires no sources . 


Following is a proof that a single source placed at any point within 
a sphere has no interference drag with the optimum distributions of 
lifting elements alone in the sphere. 


An optimum distribution of the total lift, L, within a sphere of 
radius "R" (center at the origin) is given by (see Appendix VI-l) 


^opt ~ 


rt% 2 fe 2 - r 2 



where r = Spherical radius to 
any point. Let a source loca- 
ted at an arbitrary point, P, 
within the sphere be denoted 
by S, and let P' be the pro- 
jection of P on the horizon- 
tal (x-y) plane. The poten- 
tial of S is identical with 
that caused by some lifting 

element distribution, l 

on the line between P and P 1 plus a source S' at P*. (See Ch. Ill 
Section B. The shells which have sources and sinks on the top and bot- 
tom faces respectively are arranged to form a vertical column of infini- 
tesimal cross section.) The distribution has zero net lift. 
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The interference between T-q^ and S is equal to the interference^ 

between 2 0 pt and S' plus the interference between 2 0 pt and 2^. The first 
component is zero because of the symmetry of 2 Q p^. about the x-y plane , 

(see the discussions of interference in Ch. IV). If the second compo- 
nent were not zero it would be possible to obtain a distribution of 
lifting elements alone with lower drag than 2 0 p^.. Since 2 0 p^. has the 

minimum drag by definition, the second interference component is also 
zero. This completes the proof for a particular 2 0 p^.. 

This proof can be extended to the entire family of optimum lift 
distributions in the sphere as follows. As previously mentioned, all 
of the optimum distributions produce identical effects far out on the _ 
Mach cone and far behind the wing system. Interference drag terms can 
be computed from these distant effects alone. Hence a source has the 
same interference drag with each of the optimum distributions, and this 
is zero for al 1 cases since it has been proved zero for one case . 

This proves that source distributions in a spherical volume cannot 
reduce the drag attained with any of the optimum distributions of lifting 
elements alone in that volume . 

Similar methods may be applied to ellipsoids having one principal 
axis vertical, to double Mach Cones, and to many other volumes. It is 
sufficient that the volume have a horizontal plane of symmetry, and that 
the vertical lines connecting all points in the volume with this plane 
are entirely contained within the volume. 


C . THE NON-INTERFERENCE OF SIDEFORCE ELEMENTS WITH OPTIMUM 

DISTRIBUTIONS OF LIFTING ELEMENTS IN A SPHERICAL SPACE 

As shown by Hayes there is, in general, interference between 
non-planar distributions of lifting elements, sideforce elements and 
sources . It has been proven in Ch. IX-B that there isno interference 
between a source and the optimum distribution of lifting elements alone 
in a spherical space. It remains to show a similar result for the inter- 
ference. of a sideforce element with the same optimum lift distributions. 
The proof will be carried out in a manner similar to that of the previous 
proof . 

Consider an optimum distribution of total lift L in a sphere of 
radius R; the lift distribution is given by 


^opt "" i : 

rt^R 2 - r 2 
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where r is the radius to any point from the center of the spherical 
space. Let S be a sideforce element at a point P within the sphere, 
and let P 1 be the projection of P on the xy plane. 


j? 



Fig. 9c-l 

As part of the proof it is necessary to show how a sideforce ele- 
ment can be transferred from one point to another along a line para] lei 
to the y axis. The procedure is shown in Fig. 9c- 2. First a vortex 
ring of infinitesimal, height and finite width, d, is superimposed on the 
original sideforce element; the strength and placement of the former is 


2 



jtt- 
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ELEMENT 


Fig. 9c-2 


to be such that the sideforce at P(x,y,z) is just canceled. The poten- 
tial for the vortex ring can be found by integrating the potentials for 
constant-strength infinitesimal vortex rings (Ch. IIIA.) distributed along 
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the line x, z = Constants. The second step is to superimpose on the 
vortex ring a finite width, lifting line "doublet," This latter singu- 
larity is formed by taking the limit as two equal and opposite strength 
finite width lifting elements are brought together keeping the product^ 
of lifting element strength and distance apart constant. The potential 
of the original sideforce element at P(x,y, z) plus the two added ele- 
ments, 0g(x,y, z) + is identical to that for a sideforce element 

at P(x,y-d,z). 

Thus the potential of a sideforce element "S" inside the spherical 
space is the same as that for a finite vortex ring "V, " plus a lifting 
line doublet "D, " plus a sideforce element "S'" in the vertical plane of 
symmetry (at P' in Fig. 9c-l) . The interference drag between the opti- 
mum lift distribution and- S is equal to the interference drag 

between l G pt and S' plus that between ^opt and V plus that between l opt 
and D. The last of these drags must be zero since D is a lift distri- 
bution having zero net lift; if this were not zero D could be combined 
with l C pt to form another distribution having less drag than l 0 pt 
(Section 4 h ) , in contradiction of original assumptions . Since V can be 
thought of as built up from distributions of infinitesimal vortex rings, 
which in turn are made up of source-sink doublets with axes aligned with 
the stream direction, the interference drag between V and 2 0 pt is zero 
by the proof given in Ch. IX-B. 

The only possible interference drag with Z 0 pt could be that of the 
sideforce element S' in the vertical plane of symmetry, and this can be 
shown to be zero because of the symmetry. Consider the interference 
drag of S' with lifting elements in the rear' Mach cone of S’ as shown 
in Fig. 9c- Ja. The interference drag will be due to the downwash of S’ 



s 

\ 


(a) X (6) 


Fig, 9c-3 



NACA TM lk-21 


101 


acting on the lifting elements (see Ch. IV- J); for each lift element 
which receives a downwash from S’ there is another lift element I a of 

the same strength which receives an upwash of equal magnitude . Hence 
the interference drag of S' with each pair of lifting elements in its 
rear Mach cone is zero. Consider now Fig. 9c-3b representing S’ and a 
pair of lift elements in its fore Mach cone. The interference drag 
here is due to 6idewash fields from the lift elements acting on S' . But 
for every lift element l a producing a sidewash v a on S' there is a sym- 
metrically placed 1-^ producing a sidewash v-^ = -v a on S' . Again the 

interference drag is zero. Thus there is no interference drag between S 
and Zopt and hence none between S and l 0 pt f° r this particular Zopt* 
Following the same type of reasoning as is given in Ch. IX-B, this proof 
can be extended to all optimum distributions within the sphere j this is 
so because of the uniqueness of the optimum external flow field. 

Thus it is proven that sideforce distributions , as well as source 
distributions, in a spherical space cannot reduce the drag attained with 
any of the optimum distributions of lifting elements alone in that space 

Similar methods may be applied to other spaces if those spaces have 
both a horizontal and a vertical plane of symmetry containing the free 
stream direction and meet certain convexity requirements. The latter 
can be stated as requirements that straight lines from each point within 
the space which extend to the planes of symmetry and are perpendicular 
to them must lie entirely within the space. 


Ih INTERFERENCE PROBLEMS IN CERTAIN SPACES 
BOUNDED BY MACH ENVELOPES 

Let some region "R" be chosen in the y-z plane, which is perpen- 
dicular to the flow direction. Consider the space ”S" consisting of 
points such as "P" whose fore or aft Mach cones intersect areas in the 



A set of parallel Mach planes cutting this source distribution determines 
an equivalent lineal source distribution according to the method of Hayes. 
For convenience, this equivalent lineal source distribution will be denoted 
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y-z plane which are completely contained in the region "R." An optimum 
distribution of lift in this space is given by a uniformly loaded cascade 
of airfoils (of infinitesimal chord and gap) covering the region "R, " 
since this gives constant downwash in the combined flow field. 

The resulting flow pattern is two-dimensional within the space "S." 
It then follows that the sidewash and pressure are Zero in this space 
and sideforce elements or sources introduced in "S" have no interference 
with the optimum lift distribution. 


Eh THE IH'TEHJj'ERENCE BETWEEN LIFT AND SIDEFORCE ELEMENTS 

AND AN OPTIMUM DISTRIBUTION OF VOLUME ELEMENTS 

Consider a Sears-Haack body placed on the axis of a double Mach 
cone, and place a lifting element as shown in the illustration. The 
interference drag between body 
and lifting element is composed 
of two parts, the effect of the 
body nose on the lifting element 
and the effect of the lifting 
element on the tail of the body. 


The nose of the body corre- 
sponds to a source distribution 
and produces an upwash velocity at 
the lifting element. This causes 
negative drag. The lifting element produces a positive pressure at the 
tail of the body. This also causes negative drag so the total inter- 
ference drag is negative. (This argument, of course," applies not only 
to the Sears-Haack shape but to other shapes also.) 

The total drag of the combination is equal to the drag of the 
Sears-Haack body alone plus the drag of the lifting eJLement alone plus 
the interference drag. The drag of the lifting element alone is pro- 
portional to the square of the lift it carries. However, the Interfer- 
ence drag is proportional to the first power of the lift on the element 
and to the first power of the strength- of those sources and sinks in 
the body which sire affected by interference. The lift carried by the 
element can, therefore, always be made small enough So that the drag of 
the element alone Is less (in absolute magnitude) than the interference 
drag. Thus, the total drag of the combination can be made less than 
the drag of the Sears-Haack body alone. T 

This suggests placing elements of lift’ and sideforce in a ring 
about the Sears-Haack body, and so arranged that the force on each ele- 
ment is directed radially outward from the body. This process may be 
used to construct a central body plus cylindrical shell which has zero 
drag (see Ch. IX-F) . 
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Such a system was investigated first by Ferrari^ 2 "^ and later by 
Ferri^^, and its two-dimensional analogue is the Busemann biplane. 

It, therefore, appears that the optimum distribution in space of 
volume elements alone yields minimum drag values consistent with the 
Sears-Haack values. However, the optimum distribution of volume elements 
plus lifting and sideforce elements should give zero drag for any total 
volume . 


F. RING WING AND CENTRAL BODY OF REVOLUTION HAVING ZERO DRAG 


The theoretical minimum drag value for a distribution of thickness 
elements that has no interference with lift or sideforce elements is the 
drag of a Sears-Haack body. It has been stated in Ch. IX-E that inter- 
ferences between thickness distributions and distributions of lift or 
sideforce may provide negative drag contributions which reduce the theo- 
retical minimum wave drag of a system to zero. This Section illustrates, 
for the double Mach cone volume, a central body of revolution which, 
together with a certain distribution of radial forces on a cylindrical 
shell, has zero wave drag. The method employed here to design such a 
system makes use of certain equivalences between sources and line distri- 
butions of elementary vortex shells. (These equivalences are discussed 
in Ch. III-B . ) 

Consider a radially symmetric, continuous distribution of sources 
filling a cylindrical space contained within the double Mach cone volume. 
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by F(x) . Because of the radial symmetry of this case_, F(x) is inde- 
pendent of the angle 0 on the distant control surface. A central body 
of revolution which is represented by the negative of F(x) will just 
cancel the velocities induced at the distant control surface by the 
original sources. The drag of the combined system is then zero. The 
remaining step is to relate the original source distribution to a dis- 
tribution of radial forces around the boundary of the” cylindrical space . 
It can be shown (see Ch. III-B) that a source and a sink of equal 
strength, lying on the same line parallel to the flow, direction, have 
exactly the same effect at the distant control surface as a line of con- 
stant strength elementary vortex shells connecting the source and the 
sink. If such vortex shells axe considered to replace a source distri- 
bution whose strength i6 independent of the radial distance, the forces 
on adjoining shells inside the cylinder cancel one another, while the 
forces on the outer sides of the shells next to the boundary of the 
cylindrical space determine the radial force. A cylindrical shell, haying 
this radial load distribution plus a central body of revolution which 
corresponds to the source distribution -F(x) constitute a system having 
zero drag . 

As an example, suppose that a cylinder -within a double Mach cone - 
volume is considered to contain a source distribution which varies lin- 
early with axial distance but is independent of radial distance. That 
is, the source strength per unit area inside the cylinder is 

f - -f c (i - f) 

where x is measured from the leading edge of" the cylinder, c is the 
cylinder length, and f Q is the strength of the sources at the rear face 

of the cylinder. The equivalent linear source strength corresponding 
to this original distribution is given by 
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where R is the cylinder radius. Now, if the negative of this source 
distribution is assumed to represent a body of revolution, then within 
the accuracy of slender body theory the area distribution of the central 
body is 

S(x) = - i f F(x)dx 

v J - PR 


For illustration, the dimensions of the cylinder are assumed to be 
such that pR/C = l/2; that is, the radius is half the distance between 
the axis and the apex of the double Mach cone volume. The shape of the 
central body of revolution which cancels the effect of the original 
source distribution for this case is shown in the accompanying figure. 
The distribution of radial force which can replace the original linearly 
varying source distribution is 






CENTRAL BODY OF REVOLUTION AND RADIAL 
FORCE DISTRIBUTION HAYING ZERO WAVE DRAG 
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CHAPTER X. RESULTS AND CONCLUSIONS 

It appears that certain idealized spatial distributions of lift 
and thickness may produce materially less wave drag and vortex drag than 
comparable planar systems. It is by no means certain that such advan- 
tages can be realized in practical aircraft designs, but further inves- 
tigation of specific configurations is warranted. 

One of the interesting features of spatial lift and thickness dis- 
tributions is that optimum arrangements are generally not unique. This 
may raise the problem of determining which member of an optimum family 
has the least surface area or is best adapted for structure. 

Another interesting property of spatial distributions is the inter- 
ference which may arise between lift and thickness distributions. This _ 
interference can be used to account for the zero wave drag of a Busemanh 
biplane or of Ferrari 1 s ring wing plus central body. However it is 
shown that in some cases thickness distributions have no interference 
with an optimum spatial distribution of lifting elements, and so cannot 
be used to reduce the drag due to lift in such cases. 

A number of other results are obtained in this report and detailed 
discussions of the basic singularities and Bayes 1 method of drag evalu- 
ation are included. However it is clear that the scope of the field is 
such that this investigation must be regarded as a preliminary exploration. 
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